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1.  INTRODUCTION 

The  objective  of  this  contract  is  to  obtain  data  on  the  optical 
properties  of  the  earth's  atmosphere  at  high  altitudes.  Particular  emphasis 
is  placed  on  obtaining  data  in  the  mid  infrared  region  of  the  spectrum  at  high 
spectral  resolution.  Our  group  at  the  University  of  Denver  has  a  number  of 
instruments  designed  to  be  operated  from  balloon  platforms  that  can  be  used  to 
obtain  such  data.  Data  concerning  the  optical  properties  can  be  obtained  by 
either  measuring  the  intensity  and  spectral  distribution  of  the  radiation 
emitted  by  the  atmosphere  at  high  altitude,  or  by  using  the  sun  as  an  infrared 
source  to  study  the  atmospheric  transmission  by  observing  the  spectral 
distribution  of  the  solar  radiation  at  balloon  altitudes.  The  sensitivity  of 
both  techniques  to  the  atmospheric  optical  properties  can  be  enhanced  by 
viewing  the  sun  or  the  atmosphere  close  to  the  limb.  Particular  emphasis  of 
this  program  has  been  placed  on  obtaining  atmospheric  spectral  emission  data 
from  high  altitude.  A  cryogenically  cooled  interferometer  spectrometer 
constructed  for  the  Air  Force  for  balloon  use  was  used  for  this  purpose.  This 
system,  which  has  the  acronym  SCRIBE,  has  been  flown  numerous  times  on  a 
previous  Air  Force  Program. 

These  flights  were  performed  by  the  Air  Force  balloon  detachment  and 
launched  from  Holloman  Air  Force  Base.  Reduction  of  the  number  of  groups 
using  this  capability  increased  the  per  flight  costs  and  reduced  the  number  of 
flights  which  could  be  performed  on  this  contract.  This  reduction  in  flight 
opportunities  was  partially  offset  by  an  agreement  with  NASA  to  fly  the  system 
as  part  of  the  Upper  Atmosphere  Research  Satellite  (UARS)  correlative 
measurements  program.  This  required  the  flights  be  performed  using  the 
National  Scientific  Balloon  Facility.  The  Air  Force  support  included  a  TM 
ground  station  capable  of  receiving  a  700K  bit/sec  signal.  The  NSBF  system 
did  not  have  this  capability  and  the  SCRIBE  TM  system  had  to  be  modified  so  it 
could  be  accommodated  by  the  NSBF  system.  The  other  changes  required  for  NSBF 
operation  were  minor. 

In  addition  to  the  atmospheric  measurements,  the  program  also  included 
taking  high  resolution  laboratory  data  for  compounds  of  upper  atmospheric 
interest. 


1 


2.  BALLOON  FLIGHT  PROGRAM 

The  SCRIBE  system  was  flown  on  May  23,  1989.  Since  this  flight  was 
performed  just  prior  to  the  start  of  this  contract,  reduction  of  the  data 
obtained  during  the  flight  has  been  performed  as  part  of  the  effort  on  this 
contract.  The  results  will  be  discussed  below.  Everything  worked  well  on 
this  flight  and  the  spectra  obtained  illustrated  the  capability  of  the  system. 
In  view  of  the  reduced  flight  schedule  due  to  the  increased  cost  of  flying  the 
system  it  was  decided  to  use  the  time  to  incorporate  a  second  detector  into 
the  system.  This  allows  us  to  obtain  data  from  4pm  to  8pm  in  addition  to  the 
nominal  8-15pm  region  scanned  by  the  primary  detector.  The  SCRIBE 
interferometer  is  operated  in  a  large  vacuum  dewar  so  it  can  be  cooled  to  77K 
by  liquid  nitrogen.  Incorporation  of  the  second  detector  required  extensive 
modification  to  the  vacuum  dewar  and  data  handling  system.  A  flight  with  the 
modified  system  was  scheduled  for  late  spring  1991.  The  flight  was  flown  on 
May  23,  1991. 

As  noted  in  the  Introduction,  the  additional  flight  with  the  system  was 
flown  with  balloon  support  provided  by  NASA.  The  modifications  required  to 
perform  the  flight  using  National  Scientific  Ballooning  Facility  support  were 
completed  and  the  flight  scheduled  for  February  1992  from  Barstow,  California. 
When  the  instrument  was  cooled  prior  to  shipment  to  Barstow  a  major  vacuum 
leak  occurred.  The  leak  was  associated  with  the  cryogen  fill  lines  and  only 
occurred  when  the  instrument  was  cold.  Detection  and  repair  of  the  leak 
precluded  taking  advantage  of  the  Barstow  flight  opportunity  and  the  flight 
was  delayed  until  fall  1992  from  Ft.  Sumner,  NM. 

The  final  flight  with  the  SCRIBE  instrument  was  performed  from 
Ft.  Sumner,  NM  on  October  17,  1992.  The  flight  was  launched  after 
considerable  delay  due  to  weather.  We  also  had  some  problems  with  leaks  in 
the  cryogen  fill  lines.  Our  initial  problems  with  these  lines  occurred  at  the 
fittings.  The  recent  problems  were  due  to  pin-hole  leaks  developing  in  the 
lines.  Evidently  the  extended  use  causes  some  weakening  of  the  material.  It 
is  evident  that  before  additional  flights  can  be  performed  with  the  system  the 
entire  cryogenic  system  needs  to  be  replaced.  Since  resources  were  not 
available  for  accomplishing  this,  we  decided  to  use  our  balloon-borne  cold 
grating  spectrometer  system  for  the  last  flight  on  the  program.  This  system 
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has  very  good  sensitivity  but  lower  spectral  resolution  than  SCRIBE.  The 
system  was  checked  out,  and  the  gondola  used  with  the  system  was  refurbished 
in  preparation  for  the  flight.  The  flight  with  the  system  was  scheduled  for 
the  fall  of  1994  in  support  of  a  UARS  overpass.  Preparations  for  this  flight 
were  delayed  several  days  due  to  high  surface  winds.  The  flight  was  launched 
October  9th  and  everything  went  very  well  with  excellent  data  throughout  the 
flight.  Recovery  was  accomplished  without  incident. 

3.  LABORATORY  STUDIES 

As  noted  in  the  Introduction,  a  second  technique  for  obtaining  data  on 
the  properties  of  the  upper  atmosphere  is  by  flying  infrared  solar 
spectrometer  systems.  Analysis  of  the  spectra  obtained  yield  information  on 
composition  and  physical  state  of  the  atmosphere.  During  the  period  covered 
by  this  contract  several  balloon  flights  were  performed  using  our  ultrahigh 
spectral  resolution  solar  spectrometer  system.  These  flights  were  made  using 
NASA  support,  however  the  data  obtained  are  of  interest  to  the  Air  Force. 
Phillips  Laboratory  has  a  long  history  of  developing  methods  of  predicting 
transmission  over  various  atmospheric  paths.  This  program  has  developed  a 
number  of  computer  programs  for  performing  such  calculations  including 
Lowtran,  Modtran  and  Fascode.  The  heart  of  the  calculations  with  any  of  these 
codes  is  the  molecular  line  parameters  data  base.  This  data  base,  which 
started  modestly  over  20  years  ago,  has  been  increasing  in  complexity  as  more 
spectral  regions  and  chemical  compounds  are  included  in  the  data  base. 

In  addition  to  adding  new  compounds,  the  data  base  is  continually 
updated  to  include  newer  data,  and  to  correct  mistakes  that  occurred  in  the 
earlier  versions.  This  data  base,  which  has  had  various  names  over  the  years, 
is  currently  referred  to  as  the  HiTran  Data  Base  and  the  various  versions  are 
referred  to  by  the  year  of  issue,  i.e.,  86  HiTran  Data  Base.  The  latest 
version  has  data  for  over  500,000  absorption  lines.  The  solar  spectra 
obtained  from  high  altitude  during  sunset  or  sunrise  contain  a  large  number  of 
absorption  features  which  cannot  be  observed  over  shorter  atmospheric  paths. 
The  spectra  therefore  contain  many  absorption  features  due  to  compounds 
present  at  low  concentrations  or  weak  lines  due  to  the  more  abundant 
compounds . 
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Identification  of  these  features  is  aided  by  comparison  with  the 
predictions  made  using  the  HiTran  data  base.  The  comparisons  quite  often  lead 
to  changes  in  the  data  base.  In  order  to  identify  some  features  that  do  not 
agree  with  predictions  it  has  often  proved  necessary  to  obtain  laboratory  data 
at  higher  resolution  and  over  different  temperature  and  pressure  conditions 
than  the  laboratory  data  currently  available.  During  the  period  covered  by 
this  report  we  have  obtained  such  data  for  the  following  compounds:  COCl^, 

NO^,  NO,  COF^,  HNO3,  OCS,  and  CH3CI. 

The  data  are  incorporated  into  the  HiTran  data  and  are  also  available  on 
computer  discs.  The  identification  of  the  features  in  the  spectra  are 
incorporated  into  a  series  of  Atlases  which  are  available  from  our  group. 

4.  THEORETICAL  STUDIES 

The  identification  of  weak  lines  in  the  solar  spectra  coupled  with  the 
laboratory  studies  often  enables  significant  improvements  to  be  made  in  the 
theoretical  analysis  of  many  compounds.  These  studies  play  an  important  role 
in  improving  the  HiTran  data  base.  The  results  of  these  studies  are  included 
in  the  various  versions  of  the  HiTran  data  base. 
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APPENDIX  1 


RECENT  INFRARED  EMISSION  MEASUREMENTS 
OF  PREDAWN  N2O5  OVER  NEW  MEXICO 


R.D.  Blatherwick,  F.H.  Murcray,  F.J.  Murcray,  D.G. 

Department  of  Physics,  Uni verity  of  Denver,  Denver,  CO  80208 

G.A.  Vanasse  and  M.  Hoke  nn-si 

Air  Force  Geophysics  Laboratory,  Hanscom  Field,  Bedford  MA  01731 


Introduction 

4  A  4''h  APGL  SCRIBE  system  (Stratospheric  CRyogenic 

“edaln  »205  Sn  been  completed  end  the  reeulte  ere  presented  here. 


Observations 

The  date  were  obtained  during  a  balloon  flight  >>7  the  Dni«^ 

from  Alamogordo,  DM.  Instrumentation  oonsisted  of  the  SCEXBE  LSg  cooleo 

ipi^^ptrdif& 

representing  about  5  minutes  total  observation  time. 

Analysis  and  Results 

-Tho  data  were  analyzed  by  comparison  with  synthetic  spectra  generated  using 
=™|uter Jrogram^d 

"icflaref  t%a”^m\’t:^^^ 

?K^also'rTl''orcoelfrcl^^^^^^  ro^'coJpute  the 

for  by  assuming  a  0.5%  ,  gray  background  in  lPl®‘^^fgVaround  was 

92  7“  scan,  an  additional  0.025  microwatts/cm^-sr-cm  background  was 
subLacted  from  the  data  to  achieve  agreement  with  ^he  computed  spectrum. 

The  pressure-temperature  profile  used  in  'Afloat  altitude 

Sop-sonde  releaLd  from  the  payload  shortly  after  reaching  floa^^^^ 

-0  p«|rv)v,^^rSeTH:  rriliTas%ased^n  thi 

1985  WMO  report  [5].  .u  vwm  ,.r,ioinated  (solid  line)  and 

Figures  la-lc  show  comparisons  between  the  calculate  (  927®, 
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region  arise  from  CH4,  N2O,  H2O,  and  to  a  considerably  lesser  extent,  CO2 
and  O5,  while  the  broad,  smooth  emission  feature  from  1230  to.  1265  cm  is 
attributable  to  N2O5.  The  N2O5  mixing  ratio  profile  used  in  generating  the 
synthetic  spectra  in  these  figures  was  computed  by  scaling  the  one 
dimensional  model  calculations  for  predawn  N2O5  presented  in  reference  [1] 
or  (below  25  km)  values  extrapolated  therefrom.  For  the  86.6  and  90.0 
scans,  the  scaling  factor  was  1.0.  Of  course,  one  would  expect  from 
geometrical  considerations  that  the  same  N2O5  profile  would  fit  both  of 
these  scans,  although  the  90.0®  scan  is  more  heaviliy  weighted  toward  the 
balloon  altitude.  For  these  angles,  a  six-layer  calculation  was  performed, 
and  the  N2O5  mixing  ratio  was  multiplied  by  the  same  scaling  factor  in  all 
layers.  For  the  92.7®  scan,  with  a  tangent  height  of  22.6  km,  three 
additional  layers  were  added,  and  the  mixing  ratio  profile  below  29.8  km  was 
scaled  in  three  steps  to  provide  continuity  while  achieving  a  good  fit  to 
the  observed  data.  Hence,  the  layer  from  22.6  to  24.0  km  was  scaled  by  0.5, 
the  layer  from  24.0  to  26.0  km  was  scaled  by  0.7,  the  region  from  26.0  to 
29.8  km  was  scaled  by  0.9,  and  the  mixing  ratio  above  the  balloon  was  kept 
fixed  at  the  previously  determined  values. 

Figure  2  shows  the  final  N2O5  profile  used  in  generating  all  of  the 
synthetic  spectra.  As  indicated,  the  profile  above  29.8  km  was  determined 
from  fitting  the  86.6®  and  90.0®  scans,  while  the  values  below  29.8  km  were 
determined  by  fitting  the  92.7®  scan  is  such  a  way  as  to  produce  a  smooth 
continuation  of  the  profile  above  29.8  km. 

Data  analysis  on  this  flight  is  continuing,  and  additional  results,  if 
available,  will  be  presented. 
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N2O5  VOLUME  MIXING  RATIO 


Figure  2.  N2O5  volume  mixing  ratio  profile  used  in  generating  the 
synthetic  spectra  of  Figure  1. 
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Abstract — New  line  parameters  for  two  heavy  odd  nitrogen  molecules  HNOj  in  the  V5/2V, 
region,  and  CIONO2  in  the  V4  region  are  incorporated  in  the  analysis  of  high  resolution  i.r. 
atmospheric  spectra.  The  line  parameters  are  tested  and  renormalized  vs  laboratory  spectra, 
and  then  applied  to  retrievals  from  balloon-borne  and  ground-based  solar  absorption  spectra. 


INTRODUCTION 

High  resolution  (0.002-0.003  cm“‘)  Michelson  type  interferometer  systems  are  used  at  the 
University  of  Denver  to  obtain  i.r.  solar  absorption  spectra  of  the  stratosphere  during  balloon 
flights,  as  well  as  ground-based  solar  spectra.  The  stratospheric  spectra  in  the  region  700-2200  cm"' 
cover  spectral  signatures  of  O3  (including  isotopic  species),  NO2,  HNO3,  O2,  N,,  COF2,  CIONO2, 
SFj,,  and  others.  For  a  typical  study  of  new  spectral  features  observable  in  such  spectra  see  the 
discussion  of  the  6  June  1988  flight  data  by  Goldman  et  al.‘  Ground-based  spectra  in  the  3-12 /tm 
region  are  obtained  regularly  from  Mauna  Loa  (Hawaii)  and  Denver  (CO),  in  support  of  the 
Network  for  Detection  of  Stratospheric  Change  (NDSC)  and  the  Upper  Atmospheric  Research 
satellite  (UARS). 

In  support  of  the  atmospheric  observations,  numerous  laboratory  spectra  have  been  obtained 
with  the  same  interferometer  systems.  The  laboratory  spectra  cover  mostly  the  region 
700-3000 cm"'  and  have  also  been  extended  down  to  500  or  up  to  5000cm"'.  The  molecules 
covered  include  CCI4,  CCI2F2,  CCI3F,  CF4,  CHCIF2,  CH3CI,  CHCI3,  CHjO,  CjH^,  ClONOj, 
COF2,  COCI2,  COCIF,  HCOOH,  HNO4,  H2O2,  NO,,  HNO3,  N2OJ,  and  others.  Analysis  of 
stratospheric  spectral  features,  combined  with  analysis  of  laboratory  molecular  absorption  spectra, 
improves  the  spectral  line  parameters  and  leads  to  more  complete  line  identifications  and  more 
accurate  retrieval  quantification  of  atmospheric  gases. 

In  the  present  paper,  updates  of  line  parameters  and  their  application  to  the  University  of 
Denver  spectra  will  be  presented  for  the  11/r  region  of  HNO3  and  the  780  cm"'  region  of 
CIONO2. 


RESULTS  AND  DISCUSSION 

Recent  updates  of  line  parameters  for  the  HNO3  V5/2V9  bands,  beyond  the  HITRAN  92  data 
base,^  have  allowed  improved  quantitative  analysis  of  atmospheric  spectra  in  the  1 1  /t  region.  The 
first  update  was  initiated  by  A.  G.  Maki  (unpublished)  and  included  Fermi  resonance  effects  on 
the  energy  levels  but  unperturbed  intensities.  This  provided  improved  modeling  of  the  Vj  manifolds 
in  the  860-875  cm"'  region  of  the  atmospheric  spectra,  as  reported  by  Goldman  et  al.‘  -  This  set 
has  been  used  successfully  for  most  retrievals  to  date  from  high  resolution  airborne  and 
ground-based  solar  spectra.  Subsequent  update  was  based  on  Maki  and  Wells.’  who  extended  the 
analysis  with  Fermi  resonance,  by  adding  unperturbed  intensities  to  the  published  energies  and  line 
positions.  This  update  covered  both  bands,  with  improved  agreement  in  the  Vj  manifolds  and 
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2-branch,  but  the  2V5  Q-branch  was  still  poorly  modeled.  The  more  recent  update,  based  on  the 
work  of  Perrin  et  al/  extended  the  treatment  with  both  Fermi  and  Coriolis  interactions,  including 
perturbed  intensities.  With  this  update  the  modeling  of  both  the  and  the  2V9  Q -branches  is 
significantly  improved,  and  the  2 -branches  are  now  consistent  with  most  P  and  R  manifolds.  Line 
intensity  plots  for  the  wo  latest  updates  were  showm  in  Fig.  1. 

In  order  to  apply  these  data  to  atmospheric  spectra  in  a  fashion  consistent  with  previous 
retrievals  from  these  bands  and  from  the  other  major  HNO3  bands,  the  total  band  intensities  were 
normalized  to  the  measured  intensities  of  Giver  et  al,^  divided  by  a  hot  bands  factor  of  1.30,  which 
amounts  to  multiplying  the  Perrin  et  al  lines^  (showm  in  Fig.  1)  by  1.43.  The  hot  bands  are  not 
analyzed  yet,  but  are  not  significant  in  the  spectral  intervals  used  for  stratospheric  retrievals;  they 
are  thus  not  considered  in  the  line  list  for  this  study. 

Figure  2  shows  a  typical  fitting  in  the  Vj  manifolds  with  the  latest  lines  to  data  from  the  University 
of  Denver  balloon  flight  of  6  June  1988.  The  improvements  in  the  g-branches  modeling  is 
demonstrated  in  Fig.  3,  also  with  data  from  the  6  June  1988  flight.  Comparisons  with  ground-based 
spectra  to  be  presented  elsew^here,  also  show  significant  improvements  in  the  spectral  fitting. 
Comparisons  with  retrievals  with  the  previous  tw’o  line  parameters  updates  show  consistency  of 
better  than  4%  over  the  860-870 cm”*  manifolds. 

The  first  update  of  the  initial  cross-sections  for  ClONOj  (HITRAN  1986)  was  provided  by  the 
work  of  Ballard  et  al,^  w’hich  generated  cross-sections  in  the  700-1800  cm“*  at  213-296  K.  These 
cross-sections  were  included  in  the  HITRAN  92  data  base.^ 

Recent  analysis  of  the  V4CIONO2  region  (the  most  useful  for  atmospheric  retrievals)  by  Bell  et  al,® 
provided  spectroscopic  constants  for  the  three  main  bands  in  this  region,  V4  ^^C10N02,'  V4 
^’ClONOj,  and  ^^CIONO..  Based  on  this  analysis,  line  parameters  were  generated  in 

the  V4  region,  with  the  relative  intensities  of  the  three  bands  normalized  to  0.002  cm“*  resolution 
laboratory  spectra  taken  at  the  University  of  Denver,  and  absolute  intensities  normalized  to  the 
values  of  Ballard  et  al,^  with  the  continuum  removed  (providing  consistency  of  better  than  3%  with 
previous  retrievals).  The  new  line  parameters  w^ere  then  applied  to  retrievals  in  the  780.2 cm“* 
region  from  D.U.  balloon-borne  and  ground-based  spectra  taken  at  0.002-0.006  cm”*  resolution. 
The  new  line  parameters  show^ed  improved  agreement  with  observed  data,  even  though  the  fine 
structure  is  not  fully  modeled  yet.  Intensity  plots  of  the  normalized  line  parameters  are  shown  in 
Fig.  4.  Modeling  of  the  780.2  cm”*  region  from  University  of  Denver  laboratory  data  is  showm  in 
Fig.  5.  Fitting  to  stratospheric  spectra  obtained  by  University  of  Denver  during  the  6  June  1988 
flight  is  showm  in  Fig.  6. 

The  new"  line  parameters  of  HNO3  and  CIONO^  w'ere  thus  applied  in  retrievals  of  altitude  mixing 
ratio  profiles  from  sunset  spectra  obtained  during  two  D.U.  balloon  flights  made  on  6  June  1988 
(from  Palestine,  Texas,  with  float  altitude  of  36.7  km)  and  24  July  1992  (Palestine,  Texas,  float 
altitude  of  34  km).  Sunset  scans  w'ere  obtained  during  the  6  June  1988  flight  down  to  15.60  km 
tangent  altitude,  but  only  down  to  24.30  km  tangent  altitude  during  the  24  July  1992  flight;  the 
Pinatubo  aerosol  caused  loss  of  the  source  signal  below"  this  altitude.  Both  flights  covered  the  same 
spectral  interval  of  700-1350  cm”*  with  resolution  of  -0.0025  cm”*.  The  resolution  and  the  signal 
to  noise  ratio  w^ere  slightly  belter  during  the  6  June  1988  flight.  Spectral  fitting  results  were  tested 
by  two  non-linear  spectral  least  squares  programs,  one  based  on  the  w"ork  of  Niple  et  al,^  and  one 
on  Rinsland  et  al,*^  combined  with  the  onion-peeling  method  for  the  mixing  ratio  profile.  The 
agreement  betw'een  the  two  methods  is  better  than  2%  in  any  step  of  the  onion-peeling,  assuming 
the  same  atmospheric  and  spectroscopic  parameters. 

The  retrieved  mixing  ratios  are  show'n  in  Figs.  7,  8,  which  also  include  comparisons  with  pre\ious 
results— HNO3  from  D.U.,**  ATMOS, *2  WMO  (1986),*'  and  BLISS.*^  and  CIONO:  from 

ATMOS. *^  The  results  from  the  tw’o  flights  are  included  in  the  comparisons  of  CLAES-UARS 
retrievals  and  other  correlative  measurements  and  2-D  global  model  calculations,  as  reported  by 
Mergenthaler  et  al.*'^  These  comparisons  show'  that  both  HNO3  and  CIONO2  profiles  exhibits 
significant  variabilities,  not  fully  modeled  yet,  even  for  mid  latitudes.  In  particular,  shifts  in  the 
altitudes  of  the  peak  of  the  ClONOj  profiles  such  as  seen  in  Fig.  8,  also  appear  in  the  CLAES 
results.*^  The  improved  line  parsmeters  will  provide  a  more  accurate  and  consistent  quantification 
of  these  species  from  i.r.  atmospheric  spectra. 
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Fig.  4.  Line  intensity  plots  of  the  new 

atmosohcric  spectra  for  the  v.  bands  of  ^’ClONOj  ^'ClONOj  and  v,  +  v,-  v,  of  ]«=  j'"' 

positions  are  based  on  the  spectroscopic  constants  of  Bell  et  al ,  uiih  ru^versiw^f  Denver 

4  V  IQ-’  and  JL.  =  J  =  "4.  The  relative  band  intensities  are  nonnaliaed  to  a  Un  versily  ol  uen  er 
fabomor^  spilPim.  and  absolute  intensities  normalised  to  Ballard  et  al‘  values,  with  the  continuum 

removed. 
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Fie.  5.  Spectral  least-squares  fitting  to  the  QONOj  laboratory  spectrum  with  the  new  line  parameters  in 
the  V4  ”aONOj  region.  Fine  structure  is  not  fully  accounted  for  yet. 
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Fig.  6.  Sample  retrieval  of  ClONOj  by  spectral  least-squares  fitting  to  balloon-borne  solar  spectra  obtained 
during  University  of  Denver  balloon  flight  of  6  June  I9S8. 
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Abstract.  Stratospheric  volume  mixing  ratio  profiles  of  NLO^,  CH.,  and  N2O  have  been  retrieved 
from  a  set  of  0.052  cm~^  resolution  (FWHM)  solar  occultation  spectra  recorded  at  sunrise  during  a 
balloon  flight  from  Aire  sur  lAdoun  France  (44®  N  latitude)  on  12  October  1990.  The  N2O5  results 
have  been  derived  from  measurements  of  the  integrated  absorption  by  the  1246  cm~‘  band.  Assuming 
a  total  intensity  of  4.32  x  10"^’  cm'VmolecuIe  cm"-  independent  of  temperature,  the  retrieved  N2O5 
volume  mixing  ratios  in  ppbv  (parts  per  billion  by  volume,  10"^),  interpolated  to  2  km  height  spacings, 
are  1.64  ±  0.49  at  37.5  km.  1.92  ±  0,56  at  35.5  km.  2.06  ±  0.47  at  33.5  km.  1.95  ±  0.42  at  31.5  km, 
1,60  ±0.33  at  29.5  km,  1.26  +  0.28  at  27.5  km,  and  0.85  ±0.20  at  25.5  km.  Error  bars  indicate  the 
estimated  l~cr  uncertainty  including  the  error  in  the  total  band  intensity  (±20%  has  been  assumed). 
The  retrieved  profiles  are  compared  with  previous  measurements  and  photochemical  model  results. 

Key  words.  Atmospheric  composition,  middle  atmosphere,  stratosphere.  Fourier  transform  spectros¬ 
copy.  balloon  measurements,  nitrogen  oxydes.  methane. 


1.  Introduction 

The  odd  nitrogen  family  of  molecules  (primarily  NO,  NO2,  HNO3,  CIONO2, 
HO'jNO^,  and  N2O5)  is  one  of  three  key  chemically-coupled  families  important  in 
the  catalytic  destruction  of  ozone  in  the  stratosphere  (Crutzen,  1970)  (the  others 
are  the  odd  hydrogen  and  odd  chlorine  families).  Nitrogen  pent  oxide  (N2O5)  is  a 
diumally-varying  reservoir  of  the  odd  nitrogen  family,  and  the  present  work  reports 
the  accurate  measurement  of  the  profile  of  this  molecule  at  sunrise. 

Only  a  few  measurements  of  N2O5  have  been  reported  previously.  Roscoe 
(1982)  and  Evans  (1986)  reported  tentative  observations,  based  on  low-resolution 
sunrise  measurements,  which  were  followed  by  0.015  cm~^  resolution  solar  occul¬ 
tation  measurements  of  the  743, 1246  cm"^  and  1720  cm"^  bands  by  the  ATMOS 
Fourier  transform  spectrometer  (FTS)  on  Spacelab  3  (Toon  et  ai,  1986;  Toon, 
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19S7-  Rinsland  et  ai,  1989),  nighttime  emission  measurements  of  the  1246  cm  ’ 

Id  m  W  e,al.,  1988X  ^d  morning  em.ss.on  meas¬ 
urements  of  the  same  band  with  a  balloon-borne  FTS  2  ^ 

The  present  work  reports  additional  measurements  of  the  1246  cm  NjU, 

band  ob^ed  at  sunrise  with  a  high-resolution 

solar  occultahon  mode.  The  1246  cm->  band  is  ^ 

number  of  other  atmospheric  gases,  principally  CH4  and  N,  , 

Lent  molecules  such  as  CO„  H,0.  and  COF,.  To  obtah,  accurate  cotrecnons  for 
CH.  and  N,0  interferences,  the  profiles  of  both  constituents  have  also  been 
Trieved  from  the  spectra,  and  those  results  are  also  reported. 


2.  Measurements,  Data  Analysis,  and  Results 

Table  I  summarizes  the  most  important  parameters  of  the  LP^  nrtnher 
Monitor™  Atmosphere)  instrument  and  the  balloon  flight  data  of  12  October 
1990  which  are  analvzed  here.  Instrumental  details  have  been  discussed 
PeyrJt  er  al.  (1991);  note  that  a  revised  set  of  flight  parameters  are  *iere- 

LW  the  LPMA  instrument  is  a  rapid-scanning  Bomem,  Inc.  model  DA2.0 
FTS  which  used  a  HgTdCd  detector  and  a  KCl  beamsplitter  to  cover  the 
UOO  c^'  rechon  during  this  balloon  flight.  The  spectral  signal-to-nns  "otse  of  a 
sZ^Tc^n  exceeded  100  in  the  region  of  the  1246  cm"  b^d.  About  40 

aood  spectra  were  recorded  between  solar  elevation  angles  of  -5.42  and  2  47 
frnm  thi  float  altitude  of  40 1  km.  This  altitude  was  denved  from  measurements  by 
n  onwrpre^sure  which  were  converted  to  altitude  from  a  correla¬ 

tive  pressure-temperature  profile  (discussed  below).  The  effecdve 
lution  was  about  0.052  cm"’  FWHM  (full  width  half  maximum  absorption)  _^e 
field  of  view  (FOV)  of  the  instrument  (see  Table  I)  as  determmed  y  ^  m  - 
mm  at  the  detector  and  by  the  focal  length  of  the  ou^ut 
6  X  10-^  rd,  i.e.  about  2/3  of  the  angular  diameter  of  the 

tnro  a  vertical  span  at  the  limb  of  about  2.3  km  at  a  tangent  height  of  30  km.  The 


1  Measurement  parameters  for  the  12  October  1990  balloon  flight 

Parameter 

Value 

Launch  site  . 

Float  altitude 
Instrument 

Spectral  coverage 
Spectral  resolution 
Field  of  view 

Aperture 

Scan  time 

Detector 

Beam  splitter 

Aire  sur  I'Adour,  France  (43.7'  N,  359.7  E) 

40.1  ±  0.2  km 

Bomem  DA2.01  Michelson-type  interferometer 
800-1400  cm"’  (7-12.5  pm) 

0.052  cm"'  FWHM 

6  mrad 

40  mm  diameter  (unobscured) 

50  s 

HgCdTe  (cooled  to  77  K) 

KCl 
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vertical  resolution  is  also  limited  by  the  time  betw'een  consecutive  zero  path  differ¬ 
ence  crossings  (in  our  case  2  x  50s),  which  corresponds  to  a  vertical  sampling  of 
about  2.5  km  at  the  altitude  range  of  interest. 

A  pressure-temperature  profile  calculated  by  the  National  Meteorological 
Center  (M.  Gelman,  private  communication,  1991)  for  the  time  and  location  of  the 
balloon  flight  has  been  adopted  for  the  analysis.  Refracted  ray  paths  for  each 
spectrum  were  calculated  using  the  FSCATM  program  (Gallery  et  al.,  1983)  and 
the  geometric  parameters  of  the  experiment.  The  time  of  recording  the  zero  path 
difference  signal  was  used  in  computing  the  solar  zenith  angles  (Park,  1982;  Kyle 
and  Blatherwick,  1984). 

A  nonlinear  least-squares  spectral  fitting  procedure  (e.g.,  Rinsland  el  al.,  1982) 
was  used  to  retrieve  the  CH4  and  N2O  profiles.  The  effective  instrument  line  shape 
was  first  determined  by  least-squares  fitting  isolated  CH4  lines  in  a  high  signal-to- 
noise  spectrum  produced  by  averaging  all  of  the  spectra  with  zenith  angles  less  than 
90°.  The  retrieved  instrument  function  was  held  fixed  during  the  subsequent  re¬ 
trievals.  Independent  onion-peeling  retrievals  were  run  using  each  of  the  spectral 
intervals  listed  in  Table  II.  A  sample  fit  is  shown  in  Figure  1.  The  final  profiles  have 
been  obtained  by  averaging  the  results  from  the  different  intervals.  Line  parameters 
from  the  1991  HITRAN  compilation  (L.  S.  Rothman,  unpublished  results,  1991) 
were  assumed  in  the  analysis. 

Figures  2  and  3  show  comparisons  of  the  retrieved  profiles  with  in-situ  measure¬ 
ments  obtained  at  the  same  latitude  between  1982  and  1985  (Schmidt  et  al.,  1986) 
and  2-D  photochemical  model  profiles  (M.  K.  W.  Ko,  private  communication, 
1991).  Error  bars  represent  the  total  1-a  uncertainty  of  the  present  measurements, 
about  10%  for  both  gases.  The  error  limits  were  evaluated  from  sensitivity  studies 
of  the  effects  of  uncertainties  in  each  of  the  assumed  parameters.  In  these  studies, 
retrievals  are  repeated  with  each  parameter  varied  by  the  amount  of  its  estimated 
uncertainty  while  all  others  remain  fixed  at  the  assumed  values.  The  total  errors  are 
then  calculated  from  the  RSS  (square  root  of  the  sum  of  the  squares)  of  the  indi- 


Table  11.  Spectral  intervals  used  in  the  CH4  and  N2O 
retrievals 


Interval  (cm‘‘) 

Notes 

1242.00-1247.00 

CHj  profile.  N^O  below  30  km“ 

1253.00-1255.30 

CHj  profile,  NNO  below  30  km 

1262.63-1262.96 

N^O  above  30  km 

1264.62-1264.84 

N.O  above  30  km 

1270.01-1270.25 

N^O  above  30  km 

1272.73-1273.08 

N^O  above  30  km 

1273.60-1273.85 

NnO  above  30  km 

Significant  interference  by  H^O  and  CO^  absorptions. 
The  profiles  of  both  gases  were  adjusted  as  pan  of  the  fit¬ 
ting  procedure. 
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Fie.  1.  S»pl=  .0  a  Woon-born.  .ob, 

,„-.b of 92.77-. .hich co,re.po.ded .o . min.iSlcobfri  .ol»s) 
modeling  of  the  instrument  function. 


vidual  errors  The  error  sources  and  the  average  uncertainty  in  the  retneved 
volume  mixing  ratios  resulting  from  each  are  (1)  uncertaimy  ^  *e 
spectroscopic  line  parameters  (±  4%).  (2)  uncertainty  in  the  assumed  pressure  tern 
nerature  profile  (±  2%),  (3)  uncertainty  in  the  observational  geometry  (errors  in 
pointinc  ancle  and  balloon  float  altitude)  (±5%),  (4)  uncertainty  due  to  finite 
instrum^ent  signal-to-noise  (±5%),  (5)  uncertainty  in  modeling 
line  shape  function  (±  4%).  and  (6)  uncertainty  due  to  zero  signal  level  offsets  and 

weak  channel  spectra  (±  3%).  These  last  two  sources 

cal  of  Fourier  transform  spectra.  They  are  minimal  in  the  present  expenment.  Th 
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Fig-  2-  Comparison  of  the  retrieved  LPMA  CH^  profile  with  in-situ  measurements  at  44“  N  obtained 
between  1982  and  1985  (Schmidt  et  aL,  1986)  and  a  2''D  photochemical  model  profile  computed  for 
47  N  latitude  and  the  date  of  the  LPMA  measurements  (M.  K.  W.  Ko.  private  communication,  1991). 


zero  level  can  be  checked  on  the  edges  of  the  bandpass  and  at  the  bottoms  of 
saturated  absorption  features  (H2O  and  O3  lines).  Channel  spectra,  resulting  from 
the  combined  effects  of  using  a  filter  for  isolating  the  spectral  region  of  interest  and 
a  germanium  window  for  rejecting  the  visible-UV  part  of  the  solar  spectrum,  are 
observable  but  rather  weak.  On  a  relative  basis,  the  agreement  between  the  infrared 
and  in-situ  measurements  is  very  good.  The  mean  and  the  standard  deviation  of  the 
ratio  of  the  IR  volume  mixing  ratios  to  the  in-situ  data  are  1.13  ±  0.18  for  CH4  and 
1.3410.47  for  N^O.  The  larger  standard  deviation  for  N2O  reflects  the  larger 
scatter  in  the  in-situ  measurements.  About  half  of  the  offset  between  the  CH4  data¬ 
sets  can  be  attributed  to  the  long-term  increase  in  CH4  between  the  time  of  the  in- 
situ  measurements  (1983-1985)  and  our  flight  experiment  (1990).  The  rate  of 
CH4  increase  at  the  surface  is  about  1%  per  year  (e.s.,  Khalil  and  Rasmussen, 
1990). 

Analysis  of  the  spectra  to  retrieve  the  N2O5  profile  has  followed  the  procedures 
described  previously  (Rinsland  et  al,^  1989).  In  this  method,  microwindows  spaced 
about  1  cm  ^  apart  between  1210  and  1270  cm~^,  as  given  in  Table  2  of  Rinsland  et 
ai  (1989),  were  used  to  measure  the  integrated  absorption  by  the  band.  The 
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Fio  3  Comparison  of  the  retrieved  LPMA  N,©  profile  with  in-situ  measurements  at  44  N  obtained 
between  1982  and  1985  (Schmidt  et  aL  1986)  and  a  2-D  photochemical  model  profile  computed  for 
47-  N  latitude  and  the  date  of  the  LPMA  measurements  (M.  K.  W.  Ko,  private  commumcation.  1991). 
Even  at  the  lowest  level  for  which  we  have  retrieved  the  NLO  mixing  ratio,  the  agreement  between  our 
determinations  and  the  in-situ  measurements  is  compatible  with  the  error  bars  of  each  technique. 


measurements  were  made  on  ratio  spectra  computed  by  dividing  the  individual  low 
sun  spectra  by  an  average  high  sun  spectrum  obtained  by  coadding  all  spectra  with 
zenith  angles  less  than  90'.  Corrections  for  residual  atmospheric  absorption  m  the 
microwindows  were  derived  from  simulations  calculated  with  the  retrieved  CH4 
and  N,0  profiles.  Reference  volume  mixing  ratio  profiles  (e.g.,  the  compilation  of 
Smith  (1982))  were  assumed  for  calculating  the  absorption  by  other  gases  m  this 
region.  All  of  the  species  on  the  1991  HITRAN  compilation  (L.  S.  Rothman,  un¬ 
published  results,  1991)  were  included  in  the  calculations. 

Figure  4  shows  a  sequence  of  microwindow  measurements  from  several  of  the 
spectra.  The  growth  of  the  N2O5  band  absorption  with  decreasing  tangent  height  is 
readily  apparent.  Dashed  lines  show  the  fitted  100%  transmittance  level  for  two  of 
the  spectra.  Equivalent  widths  of  the  N,0,  band  have  been  obtained  by  integrating 
the  absorption  between  1226  and  126*5  crn'l  The  results  from  spectra  with  the 
same  zero  path  difference  crossing  times  agree  within  the  uncertainties  and  have 
been  averaged.  These  times  correspond  to  the  short  period  at  zero  path  difference 
for  reverse  and  forward  scans  (one  sided  interferograms)  recorded  by  the  present 
instrument. 
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Tangent 

Heiglil 


Fio  4  Examples  of  microwindow  measurements  in  the  region  of  the  1246  cm  ‘  band  of  NLO,.  The 
values  from  each  scan,  which  have  been  corrected  for  residual  absorption  by  other  gases,  are  shown  on 
the  same  scale,  but  have  been  offset  venically  for  clarity.  Refracted  tangent  heights  of  the  spectra  are 
aiven  at  right.  Dashed  lines  show  100%  transmission  levels  determined  by  fitting  the  microwindow 
measurements  from  1210  to  1226  cm  '  and  1265  to  1270  cm  . 


The  profile  of  N;Oj  has  been  derived  from  the  measured  equivalent  widths 
using  the  onion-peeling  retrieval  method  and  an  assumed  temperature-indepen- 
dent'band  intensity  of  4.32  x  lO"''  cm-'/molecule  cm-%  This  value  is  the  same  as 
adopted  previously  (Toon  et  a!.,  1986;  Rinsland  et  al.,  1989).  Recently,  Roscoe 
(1991)  reviewed  the  published  laboratory  and  atmospheric  measurements  of  N;05 
and  recommended  an  integrated  1246  cm”’  band  intensity  of  4.01  x  10 

cm”Vmolecule  cm"-.  If  this  lower  value  is  assumed,  the  published  ATMOS  meas¬ 
urements  and  those  reported  in  this  study  need  to  be  multiplied  by  a  factor  of 
1.077. 

Figure  5  presents  a  comparison  of  our  retrieved  N2O5  profile  with  previous 
published  measurements  and  two  photochemical  model  calculations  performed  for 
conditions  corresponding  to  the  ATMOS/Spacelab  3  sunrise  measurements  at 
47'  S  on  1  May  1985.  Error  bars  on  the  LPMA  profile  increase  with  altitude,  rang¬ 
ing  from  22  to  30%.  The  uncertainty  is  dominated  by  the  absolute  uncertainty  of 
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N„Q-  Volume  Mixing  Ratio  (ppbv) 

^  D 


Fig.  5.  Comparison  of  NjOj  volume  mixing  ratio  profiles  measured  at  sunrise  by  the  LPMA  and 
ATMOS  instruments,  predawn  values  retrieved  from  spectra  recorded  by  the  SIRIS  instroment 
(Kunde  et  al.  1988).  a  post-sunrise  measurement  retrieved  from  a  spectrum  recorded  with  the 
SCRIBE  instrument  (Blatherwick  et  al.,  1989),  and  photochemical  model  calculations  generated  for 
the  conditions  of  the  ATMOS  sunrise  measurements  (1  May  1985.  at  47'  S  latitude).  Error  bars  indi¬ 
cate  the  l-siama  uncertainties  of  the  measurements.  Models  1  and  2  denote  the  model  calculations  of 

M.  Natarajan  (private  communication,  1989)  and  Allen  and  Delitsky  (1990,  Fig.  3),  respectively.  The 
latter  study  reported  profiles  derived  with  three  different  values  assumed  for  the  N^Oj  absorption 
cross  sections;  the  curve  plotted  in  this  figure  (model  2)  corresponds  to  the  profile  computed  with  the 
N^O,  cross  sections  recommended  by  DeMore  et  al.  (1987). 

20%  in  the  assumed  N705  band  intensity  at  stratospheric  temperatures.  Errors  in 
the  equivalent  width  measurements  are  also  important  at  the  higher  tangent  alti¬ 
tudes  (>35  km),  where  the  N2O5  absorption  is  very  weak.  Below  about  31  km,  the 
LPMA  profile  is  in  good  agreement  with  the  other  data.  At  higher  altitudes,  the 
ATMOS  volume  mixing  ratios  and  photochemical  model  profiles  are  lower  than 
the  present  results,  although  there  is  agreement  within  the  measurement  errors.  A 
sunrise  N,05  volume  mixing  ratio  at  30  km  has  been  calculated  by  Webster  et  al. 
(1990)  from  their  simultaneous  in-situ  measurements  of  NO,  NO2,  HNO3,  O3, 

N, 0,  pressure,  and  temperature  on  13  September  1988,  at  32” .N  latitude.  Our 
measured  value  of  1.70  ±  0.35  ppbv  at  the  same  altitude  is  significantly  lower  than 
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the  computed  value  of  2.710.4  ppbv  (Webster  er  aL  1990)  (not  plotted  in  the 
figure).  Our  measurement  results  for  CH4,  N2O,  and  N2O5  are  summarized  in 
Table  III.  For  convenience,  the  profiles  have  been  interpolated  to  2  km  altitude 
spacinss.  The  error  limits  quoted  for  all  three  gases  are  believed  to  be  conservative. 


3.  Summary  and  Conclusions 

Profiles  of  NoOj,  CH4,  and  N.O  have  been  retrieved  from  0.052  cm"’  resolution 
(FWHM)  solar  occultation  spectra  recorded  at  sunrise  near  44'  N  latitude  with  a 
balloon-borne  Michelson  interferometer.  The  measured  profiles  have  been  com¬ 
pared  with  previous  measurements  and  photochemical  model  calculations.  The 
agreement  is  generally  good. 

The  measured  spectra  contain  absorption  features  of  a  number  of  additional 
stratospheric  aases,  such  as  HNO„  '‘O,.  '‘O^O'^O,  CC1,F,,  and  CC1,F.  Work  on 
the  analysis  of  these  gases  as  well  as  spectra  recorded  dunng  other  balloon  flights  is 

in  progress. 
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Table  III.  Retrieved  volume  mixing  ratio  profiles 


Altitude 

Pressure 

Temperature 

Volume  mixing  ratio 

(km) 

(mbar) 

(K) 

CHa  (ppmv) 

NriO  (ppbv) 

N.Oj  (ppbv) 

39  5 

3.0 

246.2 

0.48  (5) 

8.0  (7) 

- 

37  5 

4.0 

241.7 

0.49  (5) 

10.3(10) 

1.64  (49) 

35.5 

33.5 

31.5 

29.5 

5.3 

7.1 

9.5 

12.8 

237,0 

231.5 

226.0 

223.0 

0.68  (7) 

0.64  (6) 

0.87  (9) 

0.81  (8) 

22(2) 

28  (3) 

66  (6) 

79  (8) 

1.92  (56) 

2.06  (47) 

1.95  (42) 

1.60  (33) 

27.5 

17.5 

220.4 

0.93  (9) 

99(10) 

1.26  (28) 

25.5 

23.8 

217.9 

1.08  (11) 

165  (17) 

0.85  (20) 

23,5 

32.4 

215.6 

1.17  (12) 

209  (21) 

ppmv  =  10"'  per  unit  volume,  ppbv  =  10'"  per  unit  volume.  Values  enclosed  in  parentheses  are 
1-a  total  uncertainties  in  units  of  the  last  quoted  digit. 
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APPENDIX  4 

The  V2  and  2p2  —  V2  bands  of  ''‘^N  ’®02:  Electron  Spin-Rotation  and 
Hyperfine  Contact  Resonances  in  the  (010)  Vibrational  State 

A.  Perrin,*  J.-M.  Flaud,*  C.  Camy-Peyret,*  a.  Goldman,!  F.  J.  Murcray,! 
R.  D.  Blatherwick,!  and  C.  P.  Rinsland! 

•  Laboratoire  de  Physique  Moleculaire  el  Applicaiions,  ^  CNRS,  Vnivershe  P.  ei  M.  Curie,  Bie  76, 

4  Place  Jussieu,  75252  Paris  Cedex  05.  France;  ^Department  of  Physics,  University  of  Denver, 
Denver,  Colorado  80208;  and  iAimospheric  Sciences  Division,  NASA  Langley  Research  Center, 

Hampton,  Virginia  2368L000} 

High-resolution  Fourier  transform  spectra  covering  the  720-920  cm"’  spectral  region  have 
been  used  to  perform  a  reanaJysis  of  the  V2  band  ((OIO)-(OOO)  vibrational  transition)  together 
with  the  first  analysis  of  the  Ivi  -  t':  hot  band  of  nitrogen  dioxide  ((020) -(010)  vibrational 
transition).  The  high-quality  spectra  show  that,  for  numerous  P:  lines,  the  hyperfine  structure  is 
easily  observable  in  the  case  of  resonances  due  to  the  hyperfine  Fermi-type  operator.  By  performing 
a  full  treatment  of  the  spin-rotation  and  of  the  hyperfine  operators,  a  new  line  list  of  the  V2  band 
(positions  and  intensities)  has  been  generated,  and  it  is  in  excellent  agreement  with  the  experimental 
spectrum.  Also,  a  thorough  analysis  of  the  2v2  —  V2  hot  band  has  been  performed  leading  to  an 
extended  set  of  new  (020)  spin-rotation  levels.  These  levels,  together  with  the  {(100),  (020), 

(001 ) }  spin-rotation  levels  deduced  previously  from  the  analysis  of  the  Pi,  2^2,  and  cold  bands 
performed  in  the  6.3-  to  l.S-pm  spectral  range  [A.  Perrin,  J.-M.  Flaud,  C  Camy-Peyret,  A.-M. 
Vasserot,  G.  Guelachvili,  A.  Goldman,  F.  J.  Murcray,  and  R.  D.  Blatherwick,  J.  Mol  Spectrosc. 

154,  391-406  (1992)]  were  least-squares  fitted,  allowing  one  to  derive  a  new  set  of  vibrational 
band  centers  and  rotational,  spin-rotation,  and  interaction  constants  for  the  { ( 100)  (020) (001 ) } 
interacting  states  of  .  ©  1993  Academic  Press,  Inc. 

I.  INTRODUCTION 

The  cold  V2  band  and  the  Ivi  —  V2  first  hot  band  of  nitrogen  dioxide,  located  at  1 3 
/iin,  correspond,  respectively,  to  the  (OIO)-(OOO)  and  (020) -(010)  vibrational  tran¬ 
sitions  of  this  molecule;  consequently,  in  order  to  have  accurate  line  positions  it  is 
necessary  to  have  precise  parameters  for  both  the  upper  and  the  lower  vibrational 
states  involved  in  these  transitions. 

For  the  (000)  and  (010)  vibrational  states,  the  most  recent  spectroscopic  parameters 
were  obtained  from  a  simultaneous  fit  of  the  available  microwave  or  double  resonance 
data  (Refs.  {1-8)  and  (9,  70)  for  the  (000)  and  (010)  state,  respectively)  and  of  the 
infixed  spin~rotation  energy  levels  obtained  from  recent  analyses  of  Fourier  transform 
spectra  recorded  in  the  far  infrared  (77)  and  in  the  13.3-mui  (72)  spectral  regions. 

In  the  analysis  of  the  (OOO)-(OOO)  band,  even- at  the  high  resolution  of  the  far- 
infrared  experimental  spectra,  the  hyperfine  structure  is  usually  not  detectable,  except 
in  cases  of  local  resonances  due  to  the  Fermi  contact  hyperfine  operator.  Accordingly, 
the  line  list  (77,  13,  14)  of  the  pure  rotation  region  of  NO2  ('^0-200  cm“^)  was 
generated  by  explicitly  taking  into  account  the  spin-rotation  and  hyperfine  structures. 

On  the  other  hand,  this  effect  could  not  be  observed  so  accurately  in  the  13.3-Mm 
experimental  spectrum  which  was  used  in  Ref.  (72)  for  the  analysis  of  the  vi  band, 
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and  consequently  the  line  list  (12-14)  generated  in  the  13.3-Mm  spectral  range  did 

not  include  the  hyperfine  structure.  were  ob- 

For  the  (020)  vibrational  state,  the  most  recent  spectroscopic  ^^^20^^^ 

tained  from  a  recent  analysis  (75)  of  the  ,  Iv.,  and  1/3  bands  (( 

( 000 )  and  ( 00 1 )  -  ( 000 )  vibrational  transitions,  respectively )  using  Founer  transfo 
pect  ;  ( 75  76)  covering  the  6.3-Mm  region.  Actually  the  (001 )  state  is  involved  in  a 
smong  Coriolis  interaction  with  the  (020)  and  (100)  states  and  Jhevery  weak  2. 

band  borrows  its  intensity  from  the  strong  band.  In  fact,  only  ^be  a  H  the  AT  = 
for  46  (because  they  appear  in  a  clear  window  between  v,  and  1/3)  and  the  Ka 
6  subband  for  19  <  <  62  (because  the  corresponding  <020 Revels  are  ^ron^ 

resonating  with  levels  of  (001 )  state)  could  be  analyzed  at  6.7 
it  must  be  pointed  out  that  for  the  energy  level  calculations  of  the  {<^0)  (020), 
(001)1  interacting  vibrational  states,  both  the  spin-rotation  and  the  vibrauonal  Co- 
riolis-type  interactions  were  explicitly  uken  into  account,  but  that  the  h>perfine  stru  - 
ture  was  not  considered: 

(i)  For  the^-type  1^3  band,  which  is  by  far  the  strongest  band  at  6.3  ^m,  the  h^erfine 
re«  affect^nly  some  .3  nonresolved  doublets,  the  spin 

in  the  upper  and  lower  levels  of  the  transition.  In  this  case,  the  bype^ine  structme  « 
blended  within  the  total  line  width,  and  the  hyperfine  component  cannot  be  obsem  . 

( ii )  Because  of  the  weakness  of  the  and  2.2  bands  the  hyperfine  structures  cannot 
be  Wisely  measured  on  the  spectra  and  consequently  the  hyperfine  resonances  were 

not  considered  for  these  bands.  ^  j  •  i  ^  -i  „m  cr.f*ctral 

Using  new  high-resolution  Fourier  transform  spectra  recorded  in  Ae  1 3-3-Mm  si^ral 

region,  we  report  in  this  paper  a  reanalysis  of  the  .2  band  (with  a  of  the 

h^erfine  structure),  together  with  a  first  analysis  of  the  hot  21-2  -  "2  band  of  N  O2. 

II.  EXPERIMENTAL  DETAILS 

The  spectra  were  recorded  in  the  720-960  cm"'  spectral  range  at  0.002  cm'*  res¬ 
olution  with  the  BOMEM  Fourier  transform  spectrometer  V^^Wfur^  of 
The  optical  path  length  was  1  m,  and  two  spectra  were  recorded  at  Pressm^  of 
1  and  2  Torr.  The  calibration  of  the  spectra  was  performed  by  means  of  N2O  absorption 
lines  (75.  19)  in  the  1150-1250  cm'*  region,  and  the  accuracy  of  the  positions  of 
unblended  lines  is  estimated  to  be  ~0.0002  cm  ' . 

III.  ANALYSIS  AND  LINE  POSITION  CALCULATIONS 

The  V2  Band 

The  analvsis  was  started  by  comparing  the  experimental  line  positions  to  the  line 
list  calculated  in  Ref.  (72).  Overall,  the  agreement  was  excellent  except  for  transiuons 
in  which  the  upper  or  the  lower  state  are  involved  in  a  hyperfine  Fermi  resonanc  . 
Wsequently,  a"new  line  list  was  generated  for  the  .2  band,  talang  into  accoum  no 
only  the  electron  spin-rotation  interaction,  but  also,  as  descnbed  m  Ref.  (77)  for  the 
(000)  state,  the  magnetic  hyperfine  hamiltonian  (i.e.,  the  Ferrni 
Uher  with  the  spin-spin  dipolar  interaction  operators)  and  the  nuclear  quadruple 
fnteraction.  For  this  calculation  we  used  the  vibrational  band  center,  rotational  spin 
Son,  and  h>parfine  consun.s  of  Rof.  (12)  and  Ref.  UD.  respecuvely,  for  .he 

unner  and  lower  State  of  the  transition.  u 

?he  agreement  this  time  proved  to  be  excellent  for  all  the  lines  However  we  obsemd 
a  -0.0004  cm”'  global  shift  of  the  observed  line  positions  relative  to  the  calcula 
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ones  which  is  due  to  a  slishtly  different  calibration  of  the  spectra  u«d  in  the  present 

worM;S.f9)compar.^.ototused(2»,in^^ 

band  center  of  the.  V2  band- £,(010)  ~  v 

°^Trs?owTe°qSity  of  thl  exlirimental  spectra  and 

in  the  Dresent  work,  we  show  in  Figs.  1  and  2  portions  around  771.8  and  7o8.0  cm 
of  tie  experimentai  spectra  together  with  line-by-line  calculations  performed  using 
ehhir  the  Une  list  generated  in  Ref.  ( 72)  ( no  hyperfine  structure )  or  the  line  list  denved 
here  (full  treatment  of  the  hyperfine  operators).  In  both  regions,  the  improvemen 

brought  by  the  new  calculation  is  very  clear.  . 

In  ttie  first  spectral  region,  the  complex  hyperfine  structure  appears  clearly  in  the 
N i  eveu  .  subband.  A.  77 1 .30  and  77, .64  cm-  i.  ,s  due  .0  a  strong  hy^rhn 
Lnance  involving  the  (0,0)[22  2  20)  and  (010)|20  2  ,8)  ™ 

cm"*  since  the  resonance  is  weaker,  one  observes  mainly  a  broadening  o 
18  doublet  structure.  The  same  hyperfine  effect  can  be  clearly  observed  in  Fig^ 
around  758  cm"'  for  lines  belonging  to  the  subband.  It  is  worth  stressing  tha 

in  both  cases  the  new  calculation,  which  involves  the  hyperfine  operators,  reproduces 
the  experimental  spectrum  with  high  accuracy. 


The  2v2  V2  Band 

Lines  from  the  21,,  -  v,  band  were  already  observed  in  the  experimental  spectnim 
( 7  7)  which  was  used  in  1 988  for  the  analysis  of  the  1 3.3-Atm  band  of  nitrogen  dioxi  ^ 

‘HievSbccausc.hbho,bu„disweu.audb^usc.hc(0^ 

available  at  that  time  were  not  precise  enou^  ( only  the  7^,  0  _  (  ) 

derived  from  the  analysis  of  the  2i,2  band  performed  at  6  6  m  (77)),  only  a  few  lines 
involving  the  Ka  =  0,  1,  2  (020)  rotational  senes  could  be  analyzed  at  13.3 

The  (020)  parameters  were  recently  significantly  improved  from  a  new  an  ysis 
the  V3}  interacting  bands  (75).  Indeed,  more  lines  involving  the  K.'OA, 

2  series  were  assigned;  also,  lines  involving  the  Tf,  -  6  rotational  senes  of  (0-  ), 
which  borrow  their  intensities  from  1,3  through  a  Conolis  interaction,  were  ob- 

toLdtidy,  from  .he  ((100),  (020).  (001 )}  and  (010)  parameter  otoiued  in 
Ref.  (75)  and  Ref.  (72)  respectively,  it  has  been  possible  to  better  predict  the  _i'2 

y,  line  Dositions,  as  performed  in  the  present  analysis. 

Because  ,he  2r,  -  r,  baud  is  rather  weak,  only  transitions  of  this  band  mvolvtng  N 
<  30  regional  quaumm  number  could  be  identified  in  the  spectrum.  For  .ranstttoM 
involving  the  =  0-2  and  Ka  =  6  (020)  energy  levels,  the  agreement  between  the 
observed  and  the  predicted  line  positions  was  excellent,  but  this  was  not  complete  y 
the  case  for  lines  involving  the  Ka  =  3-5  (020)  rotational  levels 

(UD  to  0  0030  cm"')  were  observed.  Using  these  new  data,  we  then  determine(d  the 
toWband  centers,  rotational,  spin-rotation,  and  coupling  consunts  of  the 
/(lOO)  (020),  (001)}  resonating  states.  As  m  Ref.  (7p)  the  Hamiltonian  matnx 
Explicitly  takes  into  account  both  the  Coriolis  and  the  spin-rotation 
this  least-squares  fit  calculation,  the  164  (020)  new  spin-rotation 
in  the  present  work  for  Ka  ^  6  and  N  <  30  were  combined  to  the  3230  ( 100)  (001 ) 

and  (020)  energy  levels  obtained  in  Ref  (75).  ^ 

Table  I  gives  the  list  of  vibrational  energies  and  rotational,  spin-rotation,  and  - 
pling  constants  deduced  from  the  fit,  together  with  their  estimated  uncertainty.  T 
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0.0  -1 — 
771.2 


771 .6 

cm-1 


Sj'TnTwftA  uking  into  account  the  hyperfine  sc^e^ 

and  calculation  is  excellent. 

corresponding  statistical  analysis  of  the  results  is  »™"  “  '' 

(020)  and  (001 ),  respectively. 

IV.  SYNTHETIC  SPECTRA 

This  section  presents  the  calculations'  which  have  been  perfomed  in-  order  to  gen¬ 
erate  a  precise  L  list  of  absorption  lines  of^N  "O,  in  the  13.3-Mm  region. 

^  Both  calculations  were  performed  for  a  reference  temperature  of  296  K  using  a  Z(296  K)  =  B  617.9 
partition  function  (this  value  includes  the  nitrogen  nuclear  degeneracy). 
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Fig.  2.  Experimental  spectrum  (bottom  trace)  of  nitrogen  dioxide  around  758  cm"^  and  calculated  spectra 
of  absorption  due  to  the  ’“‘N  ’^02  isotopic  species:  trace  51  and  52  are  for  the  ui  band,  respectively,  wiihout 
and  with  taking  into  account  the  hyperfine  structure;  trace  53  is  for  the  -  vz  hot  hand  (no  calculation 
of  the  hyperfine  structure).  All  spectra  have  the  same  venical  scale  but  are  displaced  for  clarity.  On  trace 
51,  we  have  marked  the  rotational  transitions  belonging  to  the  subband,  (black  dots  and  open 

circles  are  for  7  =  iV  +  ^  and  7  ==  A'  -  ^ ,  respectively).  The  agreement  between  observ’ation  and  calculation 
is  excellent. 

V2  Band 

For  the  calculation  of  the  p2  band  we  used  the  P2  transition  moment  operator  obtained 
in  Ref.  (27); 

(0oo)(O!O)^/_  =  _0  04272sc>.v  + 0.403  X  10"^{;V.v,^^} +0.1083  X  \0~'' {(f;,  iNy} .  (1) 

Also,  for  this  calculation  the  effect  of  the  hyperfine  structure  was  explicitly  taken  into 
account,  as  described  in  Ref.  (77),  through  the  standard  tensorial  formalism  using  a 
F  =  J  +  I  and  J  =  N  +  S  coupling  scheme. 

Finally,  30  620  line  intensities  were  calculated  with  an  intensity  cutoff  of  0. 1  X 
lO"*'^  cm"'  /  molecule  cm"*  and  the  following  ranges  of  quantum  numbers  and  energies; 

A'<65  7:;<13 

E'  <  2600  cm"'  £"  «  2000  cm"' 
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TABLE  1 

Vibrational  Energies  and  Rotational,  Spin-Roution,  and  Coupling  Coimnts  for  the 
{{ 100),  (020),  (001 )}  Interacting  Vibrational  States  of  N  O2 _ 


1  0  0 


0  z  0 


0  0  1 


E, 

A’ 

B’ 

c- 

At* 

a; 

s; 

HU 

K 

K 

Kn 

K 

Kn 

K 

QJ 


1319.76589 
8.0932898g 
0,431336855 
0.40928232j 
(0.2853925 
(-0.20710, 
(0.30118j 
(O.5I85 
(0.29875 
(0.337222 
(-0.2823o 
(0.220, 
(O.33O5 
(O.3782 
-0.2571 
0.10727 
{-0.553565 
(0.3250, 
(0.325e 
0.867 
-0.8439 
0.1834186 
(0.2306 
(•0.33276 
(-O.I843O3 
[  0.6005 

0.1678 
0,6322 
0.3769 
0.244 
0.29673 
-0.3568 


i  0.00012 
±  0.0000082 
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Note.  All  the  results  are  in  cm''  and  the  quoted  errors  correspond  to  one  standard  deviation.  The  consunis 
with  no  errors  were  held  fixed  during  the  fit. 


"tool) 1020J 

CCl H  0^0) 


(-0.294552i  ±  0.000088)  x  10* ‘ 

.  (O.7369  ±0.031)  X  10*5 


5‘{001){100) 

^{001)1100) 


The  total  band  intensity  was  found  to  be  5(oio)-(ooo)  -  0-542  X  10  cm  '/molecule 

^^Because  of  the  rather  high  proportion  of  so-called  AF  *  AJ  *  AJ^  “forbidden 
transitions  which  were  not  considered  in  the  1988  calculation  (I~),  the  present  cal¬ 
culated  lines  are  more  numerous  than  those  in  1988  (30  620  instead  of  8065  X  3, 
taking  into  account  the  (2  /  -I-  1 )  =  3  nitrogen  nuclear  degeneracy). 


2v2  ^  V2  Band 

For  the  hot  2v^  -  V2  band,  only  the  spin-rotation  fine  structure  was  considered  and 
its  vibrational  transition  moment  constant  was  determined  using  the  following  relation, 
valid  up  to  the  second  order  of  approximation: 

(020)(0i0)^f_  =  \/2  X  *  (2) 

The  2v^  -  V2  line  intensities  were  calculated  in  the  following  ranges: 
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TABLE  II 

Statistical  Analysis  of  the  Results  for  the  Energy  Levels  Calculation 


(100) 

(020) 

(001) 

Number  of  experimental 
spin-rotation  levels 

1187 

374 

1833 

0  <  6E  <  1x10'^ 

65.8% 

84.5% 

79.1% 

1x10'2  4  6E<  2x10'^ 

26.7% 

22.8% 

15.4% 

2x10-^  4  «E  <  7.0x10-2 

7.4% 

2.7% 

5.6% 

N  <  60  K'^  <  10 

£'  <  3200  cm”'  E"  ^  2600  cm”' 

Also,  in  order  to  be  consistent  with  the  pj  calculation,  the  intensity  cutoff  was  chosen 
to  be  0.3  X  cm”' /molecule  cm”'  at  296  K,  taking  into  account  the  3  /  +  1  =  3 
nitrogen  nuclear  degeneracy. 

Finally,  4388  lines  were  calculated  for  the  2!'2  ~  vi  band  and  the  sum  of  all 
the  intensities  was  found  to  be  5(020)-(oio)  "  0.2735  X  10  cm  ‘/molecule  cm 
at  296  K. 


V.  CONCLUSION 

Using  high-resolution  Fourier  transform  spectra,  a  new  analysis  of  the  p2  band  of 
'^'N  '^©2  involving  the  hyperfine  structure  has  been  performed  for  the  first  time  together 
with  the  first  extensive  analysis  of  the  2^2  pj  band  of  this  molecule,  leading  to  a  very 
accurate  representation  of  the  absorption  of  '"N  '^02  in  the  13.3-Mm  region.  It  must 
be  mentioned  that  lines  of  the  t'2  band  of  NO2  are  visible  in  the  University  of  Denver 
balloon  flight  spectra  of  June  1988  (22),  and  knowledge  of  their  parameters  (and  lines 
of  the  hot  band)  may  be  important  in  attempts  to  detect  the  weakly  absorbing  lines 
of  stratospheric  CIO  in  the  infrared  (22). 
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APPENDIX  5 

Stratospheric  HNO3  Measurements  From  0.002-cm“^  Resolution 
Solar  Occultation  Spectra  and  Improved  Spectroscopic 
Line  Parameters  in  the  5.8-^m  Region 

A.  Goldman,  F.  J.  Murcray,  R.  D.  Blatherwick,  J.  J.  Kosters,  and  D.  G.  Murcray 

Department  of  Physics.  University  of  Denver.  Colorado 


C.  P.  Rinsland 

Atmospheric  Sciences  Division^  NASA  Langley  Research  Center,  Hampton,  Virginia 


J.-M.  Flaud  and  C.  Camy-Peyret 

Laboratoire  de  Physique  Moleculaire  et  Applications.  CNRS.  Universite  Pierre  and  Marie  Curie,  Paris,  France 

Measurements  of  stratospheric  HNO3  presented  based  on  new  very  high  resolution  (0.002-cm"^ 
tull  width  at  half  maximum  (FWHM))  solar  occultation  spectra  recorded  at  32®N  latitude  with  a 
balloon-bome  interferometer  in  the  region  of  the  HNO3  1^2  band  at  5.8  and  improved  spectroscopic 
line  parameters  generated  for  this  HNO3  band  and  the  interfering  bands  of  O3  and  H2O.  The  line 
parameters  for  HNO3,  computed  from  the  spectroscopic  constants  reported  by  Maki,  are  a  significant 
improvementin  that  the  calculated  positions  are  more  accurate  than  previous  data  and  the  line  list  now 
includes  higher  J  and  transitions  needed  to  account  for  weaker  HNO3  absorption  features 
observable  in  the  stratosphenc  spectra.  Verifications  of  the  new  line  parameters  based  on  the 
stratosphenc  spectra  and  a  0.002-cm  resolution  HNO3  laboratory  spectrum  are  presented  along 
with  compansons  generated  with  the  parameters  from  the  1986  HITRAN  compilation.  Significant 
improvements  are  apparent  for  all  the  updated  bands.  Finally,  the  stratospheric  HNO3  profile 
retneyed  from  the  stratospheric  solar  spectra  using  a  nonlinear  least  squares  fitting  algorithm  and  the 
new  line  parameters  is  reported  arid  compared  with  previously  reported  measurements. 


1.  Introduction 

The  need  for  irriprovements  in  the  precision  and  accuracy 
of  remote  sounding  measurements  of  key  minor  constituents 
in  the  Earth’s  atmosphere  has  motivated  a  wide  variety  of 
recent  field  and  laboratory  investigations.  A  Michelson 
interferometer  capable  of  recording  broadband  spectra  in  the 
mid-infrared  at  an  unapodized  full  width  at  half  maximum 
(FWHM)  resolution  of  0.002  cm'^  has  recently  become 
operational  for  field  and  laboratory  use  at  the  University  of 
Denver  [Murcray  et  ai,  1990].  This  instrument  has  been 
flown  successfully  in  several  stratospheric  balloon  flights 
and  has  yielded  high-quality  measurements  of  a  number  of 
important  trace  constituents  [Goldman  et  aL,  1989^,  b]. 

The  present  investigation  focuses  on  the  quantitative 
measurement  of  stratospheric  HNO3  through  the  analysis  of 
the  absorption  by  its  intense  1^2  fundamental  band  in  the 
5.8-;i.m  region  of  solar  occultation  spectra  recorded  with  the 
new  liniversity  of  Denver  instrument.  Because  of  the  high 
resolution  and  high  signal-to-rms  noise  (>I00)  of  the  present 
spectra,  the  HNO3  profile  can  be  measured  accurately  up  to 
37  km,  an  altitude  seldom  reached  in  previous  investigations. 
Below  tangent  heights  of  30  km  the  HNO3  V2  band  features 
are  prominent  in  the  spectra,  but  their  absorption  is  over¬ 
lapped  by  strong  lines  from  the  +  1^2  and  P2  +  ^3 
combination  bands  of  O3,  less  numerous  but  occasionally 
strong  features  produced  by  H2O,  weak  lines  of  CH4  and 
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N2O,  and  lines  of  the  Av  =  1  sequence  of  solar  CO,  Despite 
the  complexity  of  the  absorption  this  spectral  region  has 
frequently  been  used  for  quantifying  stratospheric  HNO, 
from  high-resolution  solar  spectra  [e.g.,  Goldman  et  al., 
1980;  Coffey  et  at.,  1981,  1989;  Russell  et  al.,  1988],  As 
shown  by  the  present  analysis,  much  of  the  available  spec¬ 
troscopic  line  parameters  for  such  work  is  in  need  of 
improvement.  Line  parameters  derived  from  recent  high- 
resolution  laboratory  work  on  HNOj,  O3,  and  H2O  have 
been  used  and  are  shown  to  greatly  improve  the  quality  of 
retrievals  from  the  stratospheric  spectra.  Based  on  the  new 
spectroscopic  parameters,  a  stratospheric  HNO3  profile  is 
derived  from  the  new  solar  spectra  and  compared  with  other 
measurements  including  a  HNO3  total  column  amount  re¬ 
trieved  from  IR  solar  spectra  recorded  on  the  same  day  at 
Kitt  Peak,  about  1500  km  due  west  of  the  balloon-bome 
measurements  [Rinsland  et  al.,  1991].  To  further  verify  the 
HNO3  profile  results,  we  also  analyzed  the  balloon-bome 
spectra  to  retrieve  C62  volume  mixing  ratios  (VMRs)  and  an 
O3  profile  which  are  compared  with  other  data. 

2.  Observations 

The  stratospheric  spectral  measurements  were  recorded 
during  a  balloon  flight  conducted  from  the  National  Scientific 
Balloon  Facility  (NSBF)  in  Palestine,  Texas,  on  June  4, 
1990.  Spectra  were  recorded  both  prior  to  and  during  sunset 
from  a  float  altitude  of  37.0  ±  0.2  km,  as  estimated  from 
readings  of  on-board  pressure  transducers  converted  to 
altitude  based  on  the  1966  U.S.  Standard  Atmosphere  30'N 
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TABLE  I.  Geometry  for  the  June  4,  1990,  Balloon  Flight 
Spectra  Analyzed  in  This  Study 


Observation 
Altitude,  km 

Solar  Astronomical 
Zenith  Angle,  deg 

Refracted  Tangent 
Altitude,  km 

37.08 

90.64 

36.68 

36.98 

'  91.42 

35.02 

37.00 

92.20 

32.30 

37.05 

92.97 

28.52 

36.95 

93.73 

23.69 

36.76 

94.49 

18.14 

profile  and  a  NSBF  calibration.  The  geographic  coordinates 
of  the  balloon  gondola  at  sunset  were  3r25'N  latitude  and 
95°44'W  longitude.  Measurements  were  recorded  at  the  full 
instrument  resolution  (maximum  optical  path  diflference  of 
250  cm)  both  prior  to  sunset  and  during  sunset.  Geometric 
parameters,  the  observational  altitude,  the  solar  astronomi¬ 
cal  zenith  angle,  and  the  refracted  tangent  altitude,  are  given 
in  Table  1  for  the  scans  used  in  the  present  study.  In 
addition,  a  few  very  low  Sun  scans  (also  at  full  resolution) 
were  obtained,  but  due  to  the  large  change  in  solar  zenith 
angle  during  the  recording  time,  the  data  are  less  useful  for 
quantitative  analysis.  See  Murcray  et  al.  [1990]  for  a  detailed 
description  of  the  instrumentation  and  its  use  for  balloon- 
borne  measurements. 

3.  Spectroscopic  Line  Parameters 

An  updated  calculation  of  the  vi  band  of  HNO3  was 
generated  for  this  investigation  based  on  the  spectroscopic 
parameters  reported  by  Maki  [1988].  A  total  of  32,340  lines 
were  computed  with  a  maximum  J  and  K  value  of  74.  The 
lower  and  upper  wavenumber  limits  of  the  file  are  1650  and 
1770  cm“^  respectively.  The  calculated  line  positions  do  not 
differ  by  more  than  0.00001  cm“^  from  a  recent  independent 
calculation  by  A.  Maki  (private  communication.  1990).  An 
integrated  intensity  of  1087  cm“^  atm“*  at  296  K  corre¬ 
sponding  to  4.381  X  10“^’ cm~Vmoleculecm“~  at  296  K  has 
been  assumed  for  the  sum  of  all  1^2  band  transitions.  This 
value  has  been  derived  by  Maki  [1988]  based  on  the  inte¬ 
grated  band  intensity  measurement  of  1412  cm"*  atm"*  at 
296  K  [Giver  et  al,  1984].  Hot  bands  and  other  isotopic 
bands  are  assumed  to  produce  the  additional  absorption.  The 
room  temperature  tunable  diode  laser  intensity  measure¬ 
ments  by  May  et  al  [1987]  imply  a  total  band  intensity  0.97 
times  the  Giver  et  al  [1984]  value.  The  difference  between 
the  results  is  well  within  the  quoted  experimental  errors  of 
the  two  investigations.  For  a  discussion  of  the  comparison  of 
these  results  with  earlier  intensity  measurements,  see  May  et 
al,  [1987]. 

Maki  [1988]  noted  that  his  analysis  indicates  that  the  1/3 
band  of  HNO3  is  remarkably  free  of  perturbations.  This 
contrasts  sharply  with  the  spectrum  in  the  other  two  regions 
of  strong  absorption  by  HNO3  in  the  infrared.  At  11.3  fjLm 
the  and  21^9  bands  are  strongly  coupled  by  Fermi  reso¬ 
nance  [Maki  and  Wells,  1984];  at  7.5  jam  the  and  1/4  bands 
are  coupled  through  Fermi  and  Coriolis  resonances,  and  in 
addition,  levels  of  both  states  are  perturbed  by  interactions 
with  unobserved  levels  from  other  vibrational  states  [Perrin 
et  al,  1989].  These  interactions  make  accurate  calculations 
of  line  parameters  much  more  difficult  than  in  the  vj  band 
region.  Unfortunately,  at  present  the  agreement  between 


measurement  and  calculation  is  not  completely  satisfactory 
in  the  regions  of  the  perturbed  bands  [cf.  Goldman  et  al, 
1984.  1989^].  Since  both  May  et  al  [1987]  and  Maki  [1988] 
have  found  no  evidence  for  intensity  perturbations  in  their 
studies  of  the  V2  band,  we  have  neglected  a  Herman-Wallis- 
like  term  in  our  intensity  calculations. 

It  should  be  noted  that  the  previous  HNO3  1^2 
parameters  were  based  on  the  work  of  Maki  and  Wells 
[1980].  These  parameters  appear  on  the  1986  HITRAN  line 
parameters  compilation,  but  with  an  incorrect  intensity 
normalization  [Rothman  et  al,  1987].  The  new  work  pro¬ 
vides  a  corrected  intensity  normalization,  as  described 
above,  and  extends  the  calculations  to  higher  7  and  K.  Also, 
the  wavelength  calibration  scale  has  been  improved  due  to 
recent  heterodyne  frequency  measurements  [Wells  et  al, 
1981c,  b,  1983],  which  include  higher  rotational  levels  ob¬ 
tained  with  a  tunable  diode  laser  system.  According  to  Maki 
[1988]  the  calculated  line  positions  are  accurate  to  about 
±0.0005  cm"*  between  1680  and  1738  cm"*. 

A  Lorentz  air-broadening  coefficient  of  0.11  cm"*  atm  * 
at  295  K  and  a  temperature  dependence  of  the  Lorentz 
air-broadening  coefficient  have  been  assumed  for  all  HNO3 
lines.  These  average  values  were  reported  by  May  and 
Webster  [1989c]  based  on  their  tunable  diode  laser  labora¬ 
tory  measurements  of  HNO3  lines  in  the  7,5-/xm  spectral 
region.  The  IR  results  are  within  5%  of  calculated  linewidths 
derived  from  recent  millimeter  wave  pressure-broadening 
measurements  between  100  and  400  K  [Goyette  et  al,  1991]. 
Note,  however,  that  a  dependence  of  HNO3  pressure¬ 
broadening  coefficients  on  the  quantum  numbers  has  been 
measured  [e.g.,  Goyette  et  al,  1988]  and  is  predicted  by  the 
Anderson-Tsao-Cumutte  theory  calculations  performed  by 
Tejwani  and  Yeung  [1978].  The  new  air-broadening  param¬ 
eters  replace  the  constant  value  of  0.13  cm  "*  atm  "*  at  296  K 
and  a  7"°*^  temperature  dependence  of  the  air-broadened 
Lorentz  coefficient  which  were  assumed  for  all  HNO3  lines 
on  the  1986  HITRAN  compilation  [Rothman  et  al,  1987].  A 
constant  HNO3  self-broadened  Lorentz  coefficient  of  0.73 
cm"*  atm"*  at  296  K  [Goldman  et  al,  1980]  was  introduced 
in  the  new  linelist.  The  linelist,  which  is  in  the  1986  HITRAN 
compilation  format,  may  be  obtained  by  writing  to  the 
authors. 

To  verify  the  new  line  parameters,  two  room  temperature 
HNO3  laboratory  spectra  covering  the  1400-  to  1900-cm"* 
region  were  recorded  at  0.002-cm"*  resolution  with  a  22-cm 
path  length  absorption  cell  and  sample  pressure  of  0.4  and 
1.0  torr.  Broad  residual  H2O  background  absorptions 
formed  at  ambient  atmospheric  pressure  and  small,  narrow 
NO2  features  from  the  sample  (NO2  partial  pressure  less 
than  5%  of  the  HNO3  amount)  are  present  in  the  spectra  in 
addition  to  the  HNO3  absorption.  The  NO2  absorptions  are 
due  to  high  J,  R  branch  lines  of  the  strong  band,  centered 
at  1616,85  cm"*,  outside  the  HNO3  region.  The  background 
absorptions  have  been  partially  removed  by  dividing  the 
sample  spectra  by  a  background  (empty  cell)  spectrum. 

With  the  new  line  parameters  the  calculated  and  measured 
HNO3  lines  agree  very  well  in  position  and  intensity 
throughout  the  entire  region  of  the  1^2  band  (about  1680  to 
1740  cm"*).  However,  many  weak  lines  observed  in  the 
extreme  P  branch  of  the  laboratory  spectrum  are  not  in¬ 
cluded  in  the  synthetic  spectrum  calculations,  and  these 
lines  are  attributed  to  “hot  “  band  transitions  yet  to  be 
analyzed.  Fewer  hot  band  lines  are  encountered  in  the  R 
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Fig.  I .  Spectral  least  squares  fit  of  a  HNO3  laboratory  spectrum  in 
the  1731-cm~'  region. 


branch.  The  hot  band  lines  limit  the  accuracy  of  the  spectral 
fittings.  However,  in  narrow  intervals  with  less  interference, 
both  P  and  R  branch  features  yield  retrieved  gas  amounts 
close  (within  to  the  nominal  HNO3  amount  in  the  cell. 
An  example  of  such  a  fit  is  shown  in  Figure  1.  Our  results 
indicate  no  measurable  vibration-rotation  effect  on  the  line 
intensities  and  that  our  assumed  normalization  of  the  line 
intensities  is  correct.  These  conclusions  are  consistent  with 
the  results  reported  by  May  et  al  [1987].  Comparisons  of  the 
measured  spectra  and  simulations  also  indicate  that  the 
calculated  line  positions  at  the  1^2  b^nd  manifold  heads 
(overlapping  lines  y"  -  A/  or  /  +  \  ,  K ^  ^  ^ 

(/  +  1)  with  N.  i  integers)  are  not  fully  consistent  with  those 
outside  the  manifold  heads,  which  limits  the  agreement  that 
can  be  obtained  in  the  spectral  least  squares  fitting.  This 
effect  and  isolated  disagreements  between  the  positions  of 
measured  and  calculated  features  justify  a  new  spectro¬ 
scopic  constant  analysis  of  the  V2  band  in  the  future.  Also,  as 
expected,  fewer  hot  band  lines  are  observed  in  the  strato¬ 
spheric  spectra  because  of  the  low  temperatures. 

Figure  2  demonstrates  the  type  of  improvements  obtained 
with  the  new  HNO3  V2  band  line  parameters.  The  figure 
presents  comparisons  of  the  laboratory  data  near  1729  cm 
with  simulations  generated  with  both  the  HITRAN  1986 
lines  and  the  new  lines.  The  calculations  include  a  normal¬ 
ization  of  4.381/2.014  to  correct  the  1986  HITRAN  line 
intensities  to  the  value  of  4.381  x  10"’^  cm-Vmolecule 
cm"“  at  296  K.  The  improvement  is  particularly  important 
toward  the  wings  of  the  band,  above  1725  cm"^  and  below 
1680  cm''\  where  the  1986  HITRAN  linelist  does  not 
include  many  of  the  high  J ,  K  lines  which  show  significant 
absorption  in  the  atmospheric  spectrum.  Refinements  in  the 
line  positions  of  the  new  linelist  are  also  noticeable.  The 
improvements  are  less  dramatic  in  other  parts  of  the  band. 

The  ozone  line  positions  and  intensities  have  been  calcu¬ 
lated  from  the  results  reported  by  Malathy  Devi  et  al  [1987] 
and  Rinsland  et  al  [1990/?]  (see  Flaud  et  al  [1990a]  for  a 
listing  of  the  parameters  of  the  stronger  lines).  In  addition  to 
improving  the  positions  and  intensities  of  the  two  main 


Fig.  2.  Comparison  of  a  section  of  the  0.4  ton*  HNO3  laboratory 
spectrum  (bottom)  with  simulations  generated  using  the  1986  HI¬ 
TRAN  compilation  (top)  and  the  new  HNO3  V2  band  line  parameters 
(middle).  The  solid  triangles  beneath  the  measured  spectrum  mark 
features  that  are  missing  in  the  simulation  with  the  1986  HITRAN 
compilation  [Rothman  et  al,  1987]  but  included  in  the  simulation 
with  the  new  linelist.  Features  missing  in  both  simulations,  marked 
by  asterisks  beneath  the  measured  spectrum,  are  probably  hot  band 
lines.  The  spectra  are  offset  vertically  for  clarity. 

bands,  V2  +  (band  center  =  1726.5225  cm  )  and  vi  +  1^2 
(band  center  =  1796.2619  cm‘’),  the  new  linelist  includes 
the  two  strongest  hot  bands  2^2  +  j/3  -  ^2  and  2^2  +  vj  - 
1^2,  with  band  centers  at  1707.0034  and  1785.6456  cm  , 
respectively.  The  two  hot  bands  were  not  included  in  the 
1986  HITRAN  line  parameters  compilation  [Rothman  et  al, 
1987].  The  assumed  air-broadened  O3  line  widths  have  been 
calculated  with  an  empirical  expression  derived  from  a  fit  to 
air-broadened  line  widths  measured  from  room  temperature 
laboratory  spectra  [Flaud  et  al,  1990/?],  A  T  tempera¬ 
ture  dependence  has  been  adopted  for  the  O3  air-broadened 
widths  based  on  theoretical  calculations  performed  for  N2 
broadening  [Gamache,  1985]. 

Line  parameters  for  the  four  strongest  H2O  bands  in  the 
6-Mm  region  (excluding  HDO)  have  been  taken  from  the 
work  of  Toth  [1989,  1991].  The  bands  are  the  V2  and  the 
21^2  -  ^2  tjands  of  H2^^0,  the  1?.  bands  of  Hs^^O,  and  the  V2 
band  of  H2’*^0.  Line  parameters  from  the  1986  HITRAN 
compilation  [Rothman  et  al,  1987]  have  been  adopted  for 
the  weaker  H2O  bands  and  all  lines  of  HDO.  The  air- 
broadened  half  widths  on  the  1986  HITRAN  compilati^on 
have  been  adopted  for  all  H2O  lines  along  with  a  T  ® 
temperature  dependence  [Gamache  and  Rothman,  1988]. 
Line  parameters  for  other  gases  were  taken  from  the  1986 
HITRAN  compilation  [Rothman  et  al,  1987]. 

Figure  3  shows  two  narrow  spectral  intervals  of  the  flight 
data  dominated  by  absorption  lines  of  HNO3,  O3,  and  H2O, 
In  each  region  the  stratospheric  spectrum  is  compared  with 
two  simulations,  one  generated  with  the  new  spectroscopic 
parameters  and  the  other  with  the  1986  HITRAN  compila¬ 
tion  parameters  [Rothman  et  al,  1987].  This  comparison 
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Fig.  3.  The  1692  (upper  panel)  and  1726-cm”*  (lower  panel) 
regions  of  a  straiospheric  solar  spectrum  recorded  at  an  astronom¬ 
ical  solar  zenith  angle  of  93.73°  (tangent  height  of  23.69  km)  and 
simulations  of  the  data.  The  spectra  are  offset  vertically  for  clarity 
with  the  simulation  computed  from  the  1986  HITRAN  compilation 
[Rothman  et  aL,  1987]  at  top,  the  simulation  generated  from  the  new 
line  parameters  in  the  middle,  and  the  measured  solar  spectrum  at 
bottom  in  each  panel.  Thie  selected  intervals  show  spectral  features 
of  63,  H2O,  and  HNO3,  which  are  identified  beneath  the  lowest 
spectrum.  Note  that  two  closely  spaced  HNO3  features  near  1726.79 
cm"*  are  resolved  in  the  simulations  but  appear  as  a  single  feature 
in  the  observation.  This  is  an  example  of  the  isolated  discrepancies 
that  have  been  noted  between  the  measured  and  calculated  HNO3 
line  positions. 

demonstrates  the  improvements  that  have  been  achieved  in 
reproducing  the  atmospheric  absorption  of  the  three  gases 
with  the  new  set  of  line  parameters.  The  better  agreement  is 
in  part  the  result  of  the  corrected  normalization  of  the  HNO3 
line  intensities;  however,  other  improvements  are  readily 


apparent  (e.g.,  corrected  line  positions,  simulation  of  lines 
missing  in  the  1986  HITRAN  compilation  [Rothman  et  al., 
1987]).  Simulations  and  spectral  fittings  in  other  regions  also 
consistently  show  significant  improvements  resulting  from 
the  higher  accuracy  of  the  new  HNO3,  O3,  and  H2O  line 
parameters. 

The  new  HNO3,  O3,  and  H2O  line  parameters  noted 
above  have  been  included  in  the  1991  HITRAN  line  param¬ 
eters  compilation.  A  detailed  description  of  the  1991  HI¬ 
TRAN  data  base  is  in  preparation  (L.  S.  Rothman,  private 
communication,  1991). 

4.  Data  Analysis  and  Results 

Multiscan  onion-peeling  using  nonlinear  least  squares 
spectral  fitting  followed  by  simultaneous  multiscan  spectral 
fitting  has  been  used  to  analyze  the  balloon-borne  data. 
Details  and  examples  of  the  application  of  the  least  squares 
technique  to  derive  gas  concentration  profiles  from  IR  solar 
occultation  spectra  are  presented  in  various  papers  [cf. 
Niple,  1980;  Niple  et  al.,  1980;  Goldman  et  al.,  1980; 
Rinsland  et  ai,  1990a].  The  retrieval  algorithms  used  in  this 
study  allow  for  the  simultaneous  fitting  of  the  volume  mixing 
ratios  of  up  to  eight  gases  along  with  instrumental  parame¬ 
ters  (e.g.,  parameters  to  model  the  100%  transmission  level 
and  a  wavelength  shift  between  the  measured  and  calcula¬ 
tion  spectra)  from  the  analysis  of  each  spectral  interval. 
Because  of  the  complexity  of  the  absorption,  typically  two 
or  three  gases  have  been  fitted  in  analyzing  each  region. 
Also,  it  was  necessary  in  most  cases  to  model  the  absorption 
by  solar  CO  lines;  we  used  the  method  described  by 
Rinsland  et  al.  [1982]. 

The  zero  levels  in  the  observed  spectra  are  shifted  4-6% 
above  the  instrumental  zero  signals  with  some  frequency 
dependence.  Rather  than  fitting  the  zero  level  as  part  of  the 
spectral  least  squares  analysis,  which  can  yield  unreliable 
results  if  the  lines  are  weak,  we  used  completely  absorbed 
lines  (in  other  spectral  regions  or  scans)  to  establish  the 
corrections. 

A  pressure-temperature  profile  calculated  by  the  National 
Meteorological  Center  (NMC)  for  the  location,  date,  and 
time  of  the  balloon  flight  was  adopted  in  the  ray-tracing 
calculations  (M.  Gelman,  NMC,  private  communication, 
1990).  Gas  concentration  profiles  from  the  reference  compi¬ 
lation  of  Smith  [1982]  were  taken  as  initial  guesses  in  the 
HNO3,  O3,  and  H2O  profile  retrievals  and  in  the  calculation 
of  the  absorption  by  the  other,  less  strongly  absorbing  gases 
in  the  retrieval  intervals. 

The  spectral  intervals  studied  are  listed  in  Table  2.  In 
addition  to  selecting  intervals  to  retrieve  the  HNO3  profile, 
we  included  intervals  to  retrieve  the  CO2  VMR  and  the  O3 
profile.  The  results  from  spectral  regions  dominated  by  CO2, 
such  as  near  1933  cm"*,  indicate  that  335  ppmv  (parts  per 
million  by  volume)  CO2  is  consistent  with  the  data  (no 
detailed  profile  has  been  retrieved).  This  is  in  agreement  with 
the  value  of  337  ppmv.  anticipated  from  the  values  adopted 
previously  [Rinsland  et  al.,  1990a]  and  a  CO2  trend  of  1.42 
ppmv  yr"‘  measured  at  Mauna  Loa  [Thoning  et  al.,  1989]. 
The  ozone  profile,  obtained  from  intervals  dominated  by  O3 
lines  such  as  near  1803  cm"*,  is  very  similar  in  shape  to 
those  retrieved  from  the  May  1985  ATMOS/Spacelab  3  data 
at  28'^N  latitude  [Gunson  et  al.,  1990]  and  SAGE  II  monthly 
mean  measurements  at  32®N  in  May  and  July  1990  (J.  C. 


43 


OULDMAN  bi  Ab.:  1  KA 1  uar ncKi^^  riwm 


TABLE  2.  Spectral  Regions  Investigated  in  the  Present  Study 


spectral 

Region, 

cm~^ 

Target 

Molecules 

Primary  Interfering 
Molecules 

1686.00-1688.00 

HNO3 

O3,  solar  CO 

1690.45-1692.15 

..  HNO3 

O3,  solar  CO 

1692.90-1693.20 

HNO3 

O3 

1802.80-1804.60 

O3 

solar  CO 

1805.40-1807.40 

O3 

H2O,  NO,  solar  CO 

1933.60-1933.90 

CO^ 

1950.55-1950.95 

CO2 

O3,  COF2,  solar  CO 

1952.15-1952.50 

CO2 

solar  CO,  O3,  COF2 

1955.40-1955.70 

CO2 

1957.10-1957.50 

CO. 

O3,  COF2,  solar  CO 

Larsen,  private  communication,  1990),  but  our  results  are 
consistently  lower  by  about  12%. 

The  onion-peeling  analysis  was  begun  by  fitting  the  highest 
tangent  altitude  solar  spectrum  (37.08  km).  Absorption  by 
the  strongest  HNO3  manifolds  appears  clearly  in  that  spec¬ 
trum.  The  fit  from  each  spectral  region  was  used  to  retrieve 
individual  profile  scaling  factors  which  converted  the  target 
molecule’s  first  guess  VMR  into  an  absolute  value  in  every 
layer  with  a  lower  boundary  altitude  at  or  above  the  tangent 
height  of  the  spectrum.  VMRs  were  then  retrieved  for 
altitude  ranges  defined  by  the  tangent  heights  of  successively 
lower  Sun  spectra.  To  further  improve  the  derived  profile, 
the  simultaneous  multiscan  spectral  fitting  with  additional 
intermediate  layers  was  then  performed.  The  results  were 
adopted  as  the  final  profile. 

Figure  4  presents  a  sample  multiscan  least  squares  fit 
showing  the  1693-cm“*  region  which  is  dominated  by  medi¬ 
um-/  HNO3  V2  band,  P  branch  lines  and  weaker  absorption 
lines  of  O3  (mostly  from  the  1^2  +  v-x,  band).  Absorption  by 


uNrvERsmr  of  dekver  6/4/90  balloon  fucht 

37KM  90.64D  91.420  92.20D  92.97D  93.73D  94.49D 


WAVENUMBER  (cm"'  ) 

Fig.  4.  Simultaneous  muliiscan  spectral  least  squares  fitting  of 
the  balloon-borne  solar  spectra  in  the  1693-cm"'  region.  The 
spectra,  plotting  in  order  of  increasing  solar  zenith  angle,  are  offset 
vertically  for  clarity.  The  zero  and  100%  transmission  level  for  each 
spectra  are  shown  on  the  ordinate  axis.  The  astronomical  zenith 
angles  of  the  scans  are  listed  at  top. 


HNOjj  Volume  Mixing  Ratio  (ppbv) 

Fig.  5.  HNO3  VMR  versus  altitude  profile  retrieved  from  the 
fitting  results  shown  in  Figure  4  along  with  profiles  measured  near 
the  same  latitude.  The  ATMOS/Spacelab  3  measurements  are  from 
Russell  et  al.  [1988],  the  LIMS  data  are  from  Gille  et  aL  [1984],  the 
WMO  data  are  from  WMO  [1986],  and  BLISS  data  are  from  May 
and  Webster  [1989^]. 


solar  CO  lines  is  negligible  in  this  region  as  shown  by  high 
Sun  balloon  spectra  from  the  same  flight  and  the  ATMOS 
solar  atlas  [Farmer  and  Norton y  1989].  The  standard  devia¬ 
tion  of  the  residuals  (1.80%)  is  close  to  that  expected  from 
the  signal-to-rms  noise  ratio  of  the  spectra. 

Figure  5  presents  the  retrieved  HNO3  volume  mixing  ratio 
profile,  which  is  also  given  in  numerical  form  in  Table  3. 
Error  bars  are  the  estimated  1 -sigma  total  uncertainties, 
which  range  from  16  to  23%.  The  error  sources  and  the 
estimated  uncertainty  in  the  HNO3  VMR  resulting  from  each 
are  (1)  uncertainty  due  to  instrumental  effects  (the  finite 
signal-to-noise  ratio,  channeling,  and  zero  signal  level  off¬ 
sets)  ranging  from  10%  below  30  km  to  20%  at  37  km;  (2) 
uncertainty  in  the  assumed  tangent  pressures,  5%  at  all 
altitudes;  (3)  uncertainty  in  the  assumed  temperature  profile, 
1%  at  all  altitudes;  (4)  uncertainty  in  the  assumed  HNO3  line 
parameters,  10%  at  all  altitudes  (mostly  from  the  uncertainty 
in  the  absolute  band  intensity);  and  (5)  the  uncertainty  in  the 
simulation  of  overlapping  lines,  5%  at  all  altitudes.  The  total 


TABLE  3.  Retrieved  HNO3  Volume  Mixing  Ratio  Profile 


Altitude, 

Pressure, 

Volume  Mixing  Ratio, 

km 

mbar 

ppb,  10”^ 

19.1 

70.8 

1.88  ±  0,30 

21.0 

51.3 

2.99  ±  0.48 

22.8 

38.3 

4.22  ±  0.72 

24.8 

28.3 

6.46  ±1.16 

27.3 

19.7^  • 

5.46  ±  0.98 

29.5 

14.0 

4.99  ±  0.95 

31.4 

10.6 

2.78  ±  0.56 

32.9 

8.53 

2-75  ±  0.58 

34.3 

7.06 

1,95  ±  0.43 

35.8 

5.67 

0.92  ±  0.21 

37.3 

4.61 

0.44  ±  0.10 
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errors  have  been  computed  from  the  square  root  of  the  sum 
of  the  squares  of  the  individual  errors.  The  relative  errors  are 
largest  at  high  altitude  because  of  the  weakness  of  the  HNO3 
absorption. 

Figure  5  also  shows  HNO3  profiles  retrieved  from  data 
collected  at  about  the  same  latitude  as  the  present  measure¬ 
ments  and,  in  addition,  the  “most  likely  32°N  profile 
reported  in  World  Meteorological  Organization  (WMO) 
[1986,  Figures  10-13],  based  on  a  synthesis  of  balloon-borne 
measurements.  The  present  results  and  the  ATMOS/ 
Spacelab  3  29°N  latitude  zonal  mean  HNO3  profile  [Russell 
et  ai,  1988]  agree  very  well  except  between  25  and  30  km, 
where  the  ATMOS  VMRs  are  systematically  lower.  The 
largest  discrepancy  occurs  at  the  VMR  peak  near  25  km 
where  the  ATMOS  HNO3  mixing  ratio  is  4.8  ±  1.3  ppbv, 
significantly  lower  than  the  present  measured  value  of  6.5  ± 
1.2  ppbv  and  the  other  data  in  Figure  5.  The  Limb  Infrared 
Monitor  of  the  Stratosphere  (LIMS)  experiment  on  the 
Nimbus  7  satellite  measured  HNO3  profiles  from  the  time  of 
its  launch,  October  24,  1978,  until  May  28,  1979  [Gille  et  a/., 
1984].  The  reported  zonal  mean  LIMS  profile  for  May  1979 
at  28®N  latitude,  taken  from  Figure  16  of  Russell  et  al.  [1988], 
agrees  fairly  well  with  present  measurements.  Above  about 
34  km  the  LIMS  profile  decreases  more  slowly  with  increas¬ 
ing  altitude  than  the  present  profile  and  the  ATMOS  mea¬ 
surements  [Russell  et  aL,  1988].  As  mentioned  by  Russell  et 
al,  [1988],  the  LIMS  investigators  attributed  this  efifect  to  an 
instrumental  anomaly.  Jackman  et  al.  [1985]  noted  that  the 
LIMS  measurements  above  about  5  mbar  are  inconsistent 
with  HNO3  values  calculated  using  LIMS  H2O,  O3,  NO2, 
and  temperature  measurements.  The  HNO3  VMR  of  4.3  ± 
0.6  ppbv  at  30.6  km  retrieved  from  0.0005-cm''^  resolution 
spectra  recorded  by  the  balloon-borne  laser  in  situ  sensor 
(BLISS)  at  32®N  latitude  in  September  1988  [May  and 
Webster,  1989^]  is  very  dose  to  the  present  measurement  at 
the  same  altitude. 

A  HNO3  total  column  of  5.74  x  10^^  molecules  cm“^  has 
been  derived  by  integrating  the  measured  profile  from  18.14 
to  55.0  km.  This  value  can  be  compared  with  a  total  column 
retrieved  from  0.0053-cm~*  resolution  IR  solar  spectra  re¬ 
corded  on  the  same  day  from  Kitt  Peak  (altitude  2. 1  km, 
latitude  31,9°N,  1 1 1 .6°W),  which  is  located  about  70  km  west 
of  Tucson,  Arizona,  and  about  1500  km  west  of  the  Pales¬ 
tine,  Texas,  launch  site.  The  868.8-  to  870.0-cm  ^  region, 
which  contains  a  sequence  of  three  P  branch  manifolds  of  the 
HNO3  Vs  band,  was  analyzed  in  four  Kitt  Peak  solar  spectra 
[Rinsland  et  al,  1991].  The  line  parameters  from  the  1986 
HITRAN  compilation  [Rothman  et  al,  1987]  were  assumed 
along  with  a  correlative  pressure-temperature  profile  from 
NMC  global  satellite  and  radiosonde  measurements  (M. 
Gelman,  private  communication,  1990).  The  mean  retrieved 
Kitt  Peak  HNO3  total  column,  7.66  x  10^'  molecules  cm'^ 
is  1.33  times  the  total  column  deduced  from  the  balloon- 
borne  measurements.  It  should  be  noted  that  the  integrated 
HNO3  column  between  2  and  18  km  obtained  with  the 
reference  profile  of  Smith  [1982]  is  0.30  times  the  total. 
Therefore  after  accounting  for  the  30%  contribution  to  the 
total  column  from  layers  below  18  km,  the  Kitt  Peak  and  the 
balloon-borne  total  columns  are  consistent. 

5.  Summary  and  Conclusions 

In  the  present  work  we  have  reported  the  calculation  and 
verification  of  improved  HNO3  spectroscopic  line  parame¬ 


ters  in  the  5.8-/i.m  spectral  region  and  the  use  of  these 
parameters  along  with  improved  parameters  for  O3  and  H2O 
to  retrieve  a  stratospheric  HNO3  profile  extending  to  an 
altitude  of  37  km.  In  general,  the  retrieved  profile  is  consis¬ 
tent  with  previously  published  satellite  and  balloon-borne 
measurements  obtained  near  the  same  latitude.  Our  direct 
measurements  confirm  the  bias  in  the  high-altitude  LIMS 
measurements  noted  previously  [e.g,,  Gdle  et  al,  1984; 
Jackman  et  al,  1985].  The  excellent  agreement  between  the 
positions  and  strengths  of  HNO3  features  in  simulations  and 
in  the  very  high  resolution  stratospheric  and  laboratory 
spectra  (e.g.,  Figures  1  to  4)  confirms  the  quality  of  the  new 
line  parameters  and  their  suitability  for  use  in  the  quantifi¬ 
cation  of  HNO3  stratospheric  profiles  from  IR  spectroscopic 
measurements. 
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APPENDIX  6 


High-resolution  studies  of  atmospheric  IR  emission  spectra 

F.J.  Murcray,  F.H.  Murcray,  A.  Goldman 
R.D.  Blatherwick  and  D.G.  Murcray 
Department  of  Physics,  University  of  Denver 
Denver,  Colorado  80208 


1,  INTRODUCTION 

The  Network  for  the  Detection  of  Stratospheric  Change  (NDSC)  will  consist  of  at 
least  five  sites  located  at  different  latitudes.  Current  sites  include  a  site 
at  Lauder,  New  Zealand;  Mauna  Loa,  Hawaii;  Jungfrau joch,  Switzerland  and  two 
proposed  high  latitude  sites  in  the  Arctic  and  Antarctic. 

proposed  for  these  sites  includes  a  Fourier  Transform  Spectrometer  (FTS)  system 
to  obtain  very  high  resolution  solar  spectra  in  the  mid  infrared.  One  problem 
with  the  high  latitude  sites  is  that  the  sun  does  not  get  above  the  horizon 
during  a  significant  portion  of  the  year.  Atmospheric  emission  spectroscopy 
offers  an  alternative  way  of  obtaining  data  on  several  of  the  constituents  of 
interest.  The  technique  has  the  advantage  that  data  can  be  obtained  at  any  time 
of  day  or  night.  The  major  disadvantage  is  that  the  atmosphere  is  a  much  weaker 
source  of  radiation  than  the  sun.  Therefore  in  order  to  achieve  the  same 
spectral  resolution  at  a  given  signal-to-noise  ratio,  one  must  integrate  for  a 
much  longer  time.  The  time  required  depends  on  the  spectral  region  being  covered 
and  the  spectral  resolution  one  is  trying  to  achieve.  At  shorter  wavelengths  an 
hiqher  resolutions  the  integration  times  become  prohibitive.  At  longer 
wavelengths  it  is  possible  to  achieve  an  adequate  signal-to-noise  ratio  in  the 
spectrum  with  relatively  short  integration  times  by  sacrificing  a  modest  amount 
of  spectral  resolution.  The  trade  off  between  spectral  resolution  and 
integration  time  obviously  depends  on  the  measurement  objectives.  These  also 
dictate  the  instrumentation  required  to  make  the  measurement.  In  this  paper  we 
present  some  preliminary  results  obtained  with  two  different  systems.  These 
systems,  which  were  designed  with  different  experimental  objectives,  are  used  to 
illustrate  the  potential  of  the  technique. 


2  ♦  INSTRUMENTATION 
A.  ANTARCTIC  EMISSION  RADIOMETER 

The  first  system  was  designed  for  almost  automatic  operation  at  the  South  Pole. 
It  is  based  on  the  BOMEM  Michelson  interferometer  with  1  cm  resolution.  The 
Michelson  system  is  a  full  4  port  interferometer,  which  provides  ^ome  advantages 
in  stability.  A  mercury-cadmium-telluride  detector,  a  reference  blackbody,  a 
calibration  blackbody,  and  an  elevation  scan  mechanism  were  added  to  the  strippe 
commercial  instrument.  The  interferometer,  detector,,,  scan  mechanism,  and 
blackbodies  were  mounted  outside  the  building,  and  had  to  be  appropriately 
insulated  to  maintain  the  various  components  at  their  correct  operating 
temperatures.  Due  to  logistical  constraints,  the  detector  was  operated  with  a 
Joule-Thompson  cooler.  Gas  consumption  precluded  continuous  operation,  so  data 
were  recorded  for  2  one-half  hour  periods  each  day,  although  the  interferometer 
itself  ran  continuously. 
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Data  were  recorded  on  floppy  discs  by  a  personal  computer,  and  a  partial  backup 
copy  generated  on  the  hard  disk.  The  computer  ran  a  scheduling  program,  which 
controlled  the  Joule-Thompson  cooler  and  the  elevation  scan  mechanism.  It  also 
monitored  the  instrument  status  for  possible  errors.  Operator  maintenance  was 
required  once  per  week  to  change  the  floppy  disk,  monitor  gas  consumption,  check 
for  computer  generated  error  messages,  and  verify  that  the  instrument  optics  were 
free  of  snow. 

The  instrument  was  installed  in  the  Sky  lab  building  at  the  Amundsen-Scott  South 
Pole  station  and  began  collecting  data  on  December  5,  1989,  It  was  turned  over 
to  the  KSP  winter  over  technician  in  late  December.  The  system  operated  almost 
continuously  until  it  was  removed  in  late  January  of  1991.  One  of  the  instrument 
power  supplies  shut  down,  causing  data  loss  for  about  10  days.  This  was 
apparently  a  result  of  some  electrical  glitch  in  the  building,  but  the  supply 
restarted  after  being  powered  off  and  on.  Gas  consumption  for  the  cooler  was 
larger  than  expected,  so  the  observing  schedule  had  to  be  reduced  during  part  of 
the  year. 

The  instrument  was  designed  to  operate  in  the  500  to  1500  cm'^  region.  The 
primary  purpose  was  to  monitor  atmospheric  chemistry,  with  a  strong  secondary 
interest  in  the  radiation  budget.  The  instrument  was  set  up  to  view  two 
different  angles  (+15  and  +45  degrees  elevation).  An  initial  concern  was  with 
the  stability  and  calibration  of  an  interferometer  for  radiometric  purposes.  For 
this  reason,  the  instrument  was  calibrated  during  each  observing  sequence  by 
devoting  almost  half  of  the  time  to  viewing  the  calibration  blac]cbody. 


B.  NDSC  SOLAR  SPECTROMETER 

The  second  system  is  a  modified  Bruker  Instruments  IFS120  very  high  spectral 
resolution  interferometer,  under  test  for  the  Network  for  Detection  of 
Stratospheric  Change  (NDSC),  It  is  intended  to  observe  the  solar  spectrum  in  the 
mid  infrared,  and  is  capable  of  0.002  cm’^  unapodized  resolution. 

The  modifications  consist  primarily  of  removing  the  sample  and  detector 
compartments  from  the  standard  laboratory  instrument,  and  installing  the 
detectors  in  the  beamsplitter  compartment.  The  system  is  equipped  with  two 
detectors,  a  MCT  for  5  to  15  microns  and  an  InSb  for  2  to  5  microns.  A  solar 
tracking  system  is  usually  installed  to  bring  sunlight  into  the  instrument,  but 
for  the  emission  observations  it  was  removed,  and  skylight  was  directed  into  the 
interferometer  by  a  set  of  mirrors. 

Since  the  atmosphere  is  a  weak  source  of  radiation,  the  system  resolution  was 
reduced  to  0.06  cm"  (apodized)  for  the  emission  measurements,  allowing  a  larger 
field  aperture  and  thus  higher  throughput,  scans  were  also  coadded  for  about  20 
minutes  to  improve  the  signal-to-noise  ratio. 

Data  were  taken  from  Denver  on  a  number  of  clear  days  in  the  fall  of  1990.  The 
system  was  set  to  view  low  elevation  angles  to  enhance  the  atmospheric  path.  No 
attempt  was  made  to  calibrate  the  spectra  radiometrically. 
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3,  RESULTS  AND  Al^ALYSIS 
A.  ANTARCTIC  SPECTRA 


AS  we  have  only  had  the  data  since  January,  the  results  reported  here  are  from 
f limited  sample  and  are  preliminary.  A  check  of  the  instrument  calibrations 
shoief  a  suTdL  change  of  about  20%  in  the  fall,  then  all  calibrations  fall 

within  2%  during  the  winter  and  spring.  2%  is  roughly  the 
temperature  measurements,  so  there  was  no  measurable  change  in 

the  one  time  shift.  We  are  still  trying  to  determine  the  source  of  the  step 
change. 

A  sample  spectrum  is  shown  in  Figure  1.  This  particular  spectrum  was 
SepteLer  7,  1990  at  the  15  degree  elevation  view.  The  15  micron  CO  band  is  the 
prLominate  feature  of  the  spectrum.  It  has  an  unusual  appearance  here,  because 
of  the  strong  temperature  inversion  just  above  the  surface.  ^ 

persistent  in  the  winter  data,  and  complicates  fitting  the  CO^  band  theoretically 
because  the  standard  radiosonde  data  has  one  point  at  the 

at  the  inversion  height.  One  of  our  continuing  projects  is  to  try  to  extract 
temperature  profile  information  from  the  spectra.  The  features^ around  1050_ cm 
are^the  9.6  micron  ozone  band.  The  region  from  about  1250  cm  to  1500  cm 


predominately  water  vapor  and  methane.  Figure  2  shows  the  spectrum  near  140 
cm-"  along  with  a  theoretical  calculation  including  only  water 

required  200  microns  of  precipitable  water,  but  may  be  modified  slightly  if  a 
more  accurate  temperature  profile  can  be  obtained  from  the  CO^  region. 


V'.'kvcnuMbrr 


Figure  1.  Spectrum  taken  with  the  atmospheric  emission  radiometer 
from  the  South  Pole  on  September  7,  1990.  The  instrument  was  viewing 
15  elevation. 
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Figure  2.  Expanded  region  of  the  15  elevation  spectrum  of  Sept.  1 , 
1990.  Most  of  the  features  in  this  area  are  due  to  water  vapor. 


B.  VERY  HIGH  RESOLUTION  SPECTRA 

A  portion  of  an  emission  spectrum  obtain  with  the  NDSC  instrument  is  shown  in 
Figure  3.  This  covers  a  section  of  the  and  2v^  bands  of  HNO  .  The  strong 
features  are  weak  water  vapor  lines.  A  calculated  emission  spectrum  is  also 
shown,  using  a  standard  temperature  and  water  vapor  profile.  The  ordinate  scale 
applies  to  the  calculated  spectrum,  the  observed  spectrum  was  linearly  scaled  and 
a  zero  subtracted  to  match  the  calculation.  The  agreement  is  very  good  except 
for  the  line  at  87  5.7  cm‘^,  which  may  have  an  incorrect,  strength  listed  in  the 
data  base.  Another  discrepancy  is  at  876.6  cm"^  which  is  probably  due  to  HNO^. 

Figure  4  is  an  expanded  view  of  the  865  to  870  cm*^  spectral  interval.  The  upper 
spectrum  (A)  in  this  figure  is  a  solar  absorption  spectrum  taken  at  the  same 
resolution  as  the  emission  data  and  plotted  upside-down.  The  second  spectrum  (B) 
is  the  atmospheric  emission  in  the  same  interval.  The  bottom  (C)  spectrum  is  a 
calculation.  The  emission  features  are  mainly  due  to  HNO^.  The  broad  feature 
in  the  experimental  data  in  the  867  to  868  cm*^  and  869  to  870  cm*^  area  is  a 
channel  spectrum  due  to  a  non-wedged  filter.  Some  weak  HNO^  lines  appear  in  the 
observed  spectra  but  are  not  in  the  calculation.  The  emission  spectrum  also 
shows  a  feature  near  868  cm“^,  which  is  probably  due  to  NH^.  NH^  is  highly 
variable  in  concentration  near  the  surface,  and  shows  considerable  enhancement 
in  the  emission  spectrum  compared  to  HNO^  because  it  is  warmer  since  it  is  near 
the  surface. 
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Figure  3.  Portion  of  an  atmospheric  emission  spectrum  obtained  at  an 
elevation  angle  of  15*.  This  is  compared  with  a  spectrum  calculated 
using  a  standard  atmospheric  temperature  profile.  The  observed 
spectra  is  shown  on  the  dashed  line  and  has  been  adjusted  to  fit  the 
water  vapor  features  in  the  calculated  spectra. 


Wavenymber/cm-l 

Figure  4.  Enlarged  portion  of  the  emission  spectrum  given  in  Figure 
3  (Curve  B).  This  spectrum  is  compared  with  a  solar  spectrum  (Curve 
A)  obtained  with  the  same  instrument  at  about  the  same  elevation 
angle.  Also  shown  (Curve  C)  is  a  theoretical  spectrum  for  the  same 
spectral  interval. 


4.  CONCLUSIONS 


l'7e  have  presented  atmospheric  emission  spectra  obtained  with  two  different 
spectrometer  systems.  The  first  system  was  designed  for  emission  work.  Spectra 
were  obtained  under  adverse  conditions  in  the  Antarctic,  and  are  still  of  good 
absolute  accuracy.  The  second  system  demonstrates  the  sensitivity  that  can  be 
achieved  even  at  higher  spectral  resolution.  Although  the  system  was  not 

calibrated,  and  is  not  particularly  well  designed  for  emission  work,  it  still 
shows  that  mid  infrared  atmospheric  emission  spectra  can  be  obtained  with  good 
signal-to-noise  in  a  reasonable  length  of  time  at  relatively  high  resolution. 
A  properly  designed  high  resolution  system  should  achieve  high  accuracy, 
sensitivity  and  resolution,  thereby  permitting  measurements  of  many  atmospheric 
constituents  when  solar  spectra  can  not  be  obtained. 
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APPENDIX  7 

The  Fundamental  Quadrupole  Band  of  ’“Na:  Line  Positions  from  High-Resolution 
Stratospheric  Solar  Absorption  Spectra 

The  purpose  of  this  note  is  to  report  accurate  measurements  of  the  positions  of  O-  and  S'-branch  hnes  of 
the  ( 1-0)  vibration-rotation  quadrupole  band  of  molecular  nitrogen  ('“N:)  and  improved  Dunham  cwffiaents 
derived  from  a  simultaneous  least-squares  analysis  of  these  measurements  and  selected  infrared  and  far 
infrared  data  taken  from  the  literature.  The  new  measurements  have  been  denved  from  stratosphew  solar 
occultation  spectra  recorded  with  Fourier  transform  spectrometer  (FTS)  instruments  operated  at  unapodized 

spectral  resolutions  of  0.002  and  0.01  cm  r  •  te> 

The  motivation  for  the  present  investigation  is  the  need  for  improved  N2  line  parameters  for  in  IR 

atmospheric  remote  sensing  investigations.  The  5  branch  of  the  N,  (1-0)  quadrupole  band  is  id^  for 
calibrating  the  line-of-sight  airmasses  of  atmospheric  spectra  since  the  strongwt ^ '  Snmical 
atmospheric  window,  their  absorption  is  relatively  insensitive  to  temperature  and  is  moderately  strong  (typ 
line  center  depths  of  10  to  50%  in  high-resolution  ground-based  solar  spectra  and  in  lower  stratosphenc 
solar  occultation  spectra),  and  the  volume  mixing  ratio  of  nitrogen  is  constant  in  the  atmophere  and  wel 
known.  However,  a  recent  investigation  has  shown  the  need  to  improve  the  accuracnes  of  the  N.  line  turnons 
intensities,  air-broadened  half-widths,  and  their  temperature  dependences  to  fully  exploit  this  i^ibration 
capability  (/).  The  present  investigation  addresses  the  problem  of  improving  the  accuracy  of  the  Nj  h  e 

‘^SrSen  is  a  homonuclear  diatomic  molecule,  and  therefore  electric  dipole  pure  rotation  and  vibration-- 
rotation  transitions  within  the 's;  ground  electronic  state  are  forbidden.  For  this  reason, 
and  rotational  constants  of  the  ground  electronic  state  have  remained  poorly  determined  with  niost  of  the 
reported  values  derived  from  medium-resolution  Raman  spectra  recorded  in  the  far  infrared  (2  4)  and 
infrared  (i-5)  regions.  Reuter  et  al.  (6)  used  a  White  cell  to  achieve  an  absorption  length  of  434  m 
and  a  FTS  to  obtain  measurements  of  line  positions  in  the  O  and  ^branches  of  the  (  )  in  ra  lon-ro  i 

quadrupole  band.  The  positions  derived  in  Ref.  (6)  were  significantly  more  precise  than  achieved  previously 
and  the  use  of  standard  IR  calibration  techniques  yielded  a  wavenumber  calibration  supenor  to  those  oft 
Raman  investigations.  Measurements  with  line  position  accuracies  of  0.003  to  O.OlO  cm  were  rei^  e 
ta  ..  0.0um-'  r^olutta.  Th.  «« “ 

weakness  of  the  N.  lines  in  the  measured  spectrum  (maximum  line  center  absorption  of  a  few  tenths  of  one 

'^Thrlratospheric  measurements  reported  here  were  obtained  from  solar  spectra  recorded  dunrig  flights 
bv  two  different  instruments.  The  first  set  of  data  was  acquired  from  orbit  on  April  30-May  1,  1985,  W 
atmospheric  trace  molecule  spearoscopy  (ATMOS)  instrument,  a  high  speed  FTS  o^raung  at  an  unapor^^ 
spectr^  resolution  of  0.01  cm"  during  its  first  Shuttle  mission  onboard  Spacelab  3- For  an 
ATMOS  instrument  and  experiment  and  the  Spacelab  3  results,  see  the  review  paper  by  Farmer  (  7).  ATMOV 
related  data  processing  and  science  analysis  methods  are  described  in  a  paper  by  Norton  and  I^nslaiid  (8). 
The  other  set  of  measurements  was  obtained  with  a  FTS  achieving  an  unapodized  resolution  of  0.002  crn 
onboard  a  balloon  platform.  The  instrument  is  a  modified  commercial  interferometer  operated  by  the 
atmospheric  spectroscopy  group  at  the  University  of  Denver.  Solar  obse^ations  were 
to  and  during  sunset  from  a  float  altitude  of  36.6  ±  0.5  km  dunng  a  flight  from  the  National  Scientific 
Balloon  Facilitv  in  Palestine,  Texas,  on  June  4.  1990.  A  description  of  the  instrumentation  and  examples 
of  spectra  from  several  balloon  flights  can  be  found  in  the  paper  by  Murcray  et  al.{9). 

The  positions  of  unblended  N:  lines  between  6'(4)  and  S(20)  (2373  to  24^2  cm  )  have  l^n 

measured  from  ATMOS  spectra  recorded  with  a  filter  covering  approximately  1580  to  3400  cm  .  The 

continuum  signal-to-rms  noise  in  the  Nj  5-branch  region  is  -100  iri  t^tte 

N,  5-branch  lines  first  appear  above  the  noise  level  at  a  tangent  altitude  of  ~40  km  (10).  Because  the 
atmospheric  density  increases  exponentially  with  decreasing  altitude,  the  nitrogen  lines  strengthen  rapidly 
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Fig.  1.  Examples  of ‘'‘N2  (1-0)  band  quadrupole  lines  in  the  ATMOS  and  University  of  Denver  infrared 
stratospheric  solar  occultation  spectra.  The  upper  panel  (A)  shows  the  region  of  the  5(9)  line  in  a  zonal 
average  .ATMOS  spectrum  with  tangent  altitude  of  1 9.7  km.  The  slight  decrease  in  signal  to  lower  wavenumbers 
is  caused  by  continuous  absorption  produced  by  the  N2  pressure-induced  ( I-O)  band  and  sub~Lorcntzian 
absorption  by  the  overlapping  far  wings  of  the  intense  CO2  fundamental  band  (77).  The  region  of  the 
D(  16)  line  is  plotted  in  the  lower  panel  (B).  The  spectrum  was  recorded  at  a  tangent  altitude  of  23.4  km 
w'ith  the  University  of  Denver  FTS.  The  nitrogen  lines  are  identified  and  marked  with  vertical  ticks  in  both 
panels. 


with  decreasing  tangent  altitude.  The  most  accurate  measurements  are  obtained  from  spectra  recorded  near 
a  tangent  altitude  of  20  km  w^here  the  stronger  5-branch  lines  are  prominent  features  in  the  ATMOS  spectra 
(70).’  At  that  altitude,  the  atmospheric  pressure  is  ^40  Torn  so  that  the  effect  of  pressure-induced  line 
shifts  should  be  small  At  lower  tangent  altitudes,  the  signal-to-noise  ratio  is  decreased  by  continuum  absorption 
resulting  from  the  collision-induced  fundamental  vibration-rotation  band  of  N2  and  sub-Lorentzian  absorption 
produced  by  the  superposition  of  the  far  wings  of  the  intense  CO2  ^^3  fundamental  band  (77).  Also,  the 
absorption  by  other  telluric  gases  generally  increases  with  decreasing  tangent  altitude  so  that  interference 
effects  become  more  of  a  limitation.  For  the  present  w'ork.  to  improve  the  signal-to-noise  ratio,  we  derived 
the  line  positions  from  zonal  average  ATMOS  spectra  obtained  by  coadding  nearly  equal  tangent  altitude 
spectra  recorded  during  four  sunset  occultations  near  30®N  latitude.  The  line  positions  were  measured  with 
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a  second  derivative  line  finder  algorithm  and  calibrated  using  accurate  line  positions  reported  for  the  N2O 

(72)  or  j'j  +  2v2  bands  ( 7i).  Nitrogen  line  positions  from  1  to  4  ATMOS  zonal  spectra  have  been  averaged 
for  use  in  the  subsequent  analysis;  the  corresponding  uncertainties  were  estimated  on  the  basis  of  the  standard 
deviation  of  the  individual  measurements  and  the  calculated  uncertainty  derived  from  the  formula  given 
by  Reuter  ei  al.  {6).  A  minimum  uncertainty  limit  of  0.(X)02  cm"'  was  adopted  to  account  for  the  uncertainty 
in  the  absolute  wavenumber  calibration. 

The  balloon-borne  solar  spectra  were  recorded  with  a  filter  covering  the  1580  to  2230  cm"'  region  and 
hence  provide  measurements  of  only  high-y  N2  lines  in  the  O  branch.  A  total  of  6  N2  lines  between  0(13) 
and  0(19)  were  selected  for  measurement  based  on  comparisons  of  the  solar  spectra  with  synthetic  spectrum 
calculations.  All  of  the  measured  transitions  are  calculated  to  be  free  of  significant  blending  by  other  telluric 
lines.  The  0(20)  N2  line  also  appears  to  be  present,  but  it  is  at  the  limit  of  detection  and  was  not  included 
in  the  analysis.  As  the  N2  lines  occur  near  the  edge  of  the  band-pass  filter,  the  continuum  signal-to-rms  noise 
ratio  of  a  single  spectrum  decreases  from  ^100  at  the  0(  19)  line  to  '-'50  at  the  0(  13)  line.  The  positions 
of  the  N2  lines  were  measured  from  one  to  three  spectra  recorded  between  tangent  altitudes  of  28.1  and 
18.0  km.  Mean  values  were  used  in  the  subsequent  analysis.  Calibration  of  the  individual  solar  spectra  was 
achieved  by  comparing  measured  positions  of  strong,  isolated  lines  (e.g..  H2O  and  CO2  transitions)  with 
values  derived  from  stratospheric  spectra  recorded  with  the  same  instrument  during  an  earlier  balloon  flight. 
The  calibration  of  this  earlier  dataset,  which  covers  1400  to  I960  cm"' ,  is  based  on  H2O  standard  positions 
(74),  scaled  by  the  multiplicative  factor  recommended  by  Brown  and  Toth  (75).  The  frequency  scale  of 
the  new  flight  data  was  also  checked  using  accurate  NO2  band  line  positions  (76),  as  well  as  direct 
comparisons  of  measured  H2O  line  center  positions  with  the  data  of  Brown  and  Toth  (75).  A  second 
derivative  line  finder  was  used  here  to  determine  the  location  of  the  line  centers.  Uncertainties  were  estimated 
as  described  earlier.  Figure  1  presents  examples  of  the  N2  lines  in  the  ATMOS  (A)  and  University  of  Denver 
(B)  stratospheric  solar  spectra. 

The  line  position  measurements  were  analyzed  with  the  Dunham  formula  for  the  upper  and  lower  state 
term  values  (77) 

Z  5)'[W+  1)]^  (1) 

ij 


TABLE  I 

Dunham  Coefficients^  for  "‘N2 


Dunham 

This  Study 

Coefficient 

Case  A 

Case  B 

Ref,  6 

Ref.  5 

1.9982673(43) 

1.9982369(56) 

1.998286(13) 

Yo2  X  10® 

-5.761(10) 

-5.747(14) 

-5.763(40) 

Yo,  X  10=2 

0. 

4.78^ 

0. 

Yio 

2329.91170(20) 

2329.91156(21) 

2329.91239(70) 

Y.j  X  10^ 

-1.73717(9) 

-1.73030(19) 

-1.73121(67) 

1.72978 

X 

o 

0. 

-1.01(56) 

0. 

-0.832 

Yji  X  10= 

0. 

-3.28(20)' 

-3.28(20)' 

-4.15 

Weighted 

Error 

2.72  X  10"* 

2.67  X  lO"'- 

1.99  X  10“^ 

a)  Values  in  parentheses  are  standard  deviations  in  units  of  the  last  quoted 
digit  of  the  coefficient. 

b)  Calculated  from  case  A  Dunhaia  coefficients  and  Eq.  2  of  Reuter  et  al.  (6); 
value  was  held  fixed  in  the  analysis. 

c)  Taken  from  Lofthus  and  Krupenie  (1^);  value  was  held  fixed  in  the  analysis. 
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with  the  coefficients  obtained  from  a  weighted  linear  least-squares  fit  of  the  stratospheric  measurements 
and  additional  data  selected  from  reported  laboratory  experiments.  Each  weight  was  given  by  1  where 
a  is  the  estimated  uncertainty.  The  additional  data  were  required  because  the  stratospheric  dataset  covers 
only  a  limited  range  of  transitions  in  the  O  and  S  branches.  Three  additional  datasets  were  included:  the 
ground  state  Raman  measurements  of  Bendtsen  (i),  the  Raman  (1-0)  band  (9-branch  measurements  of 
Bendtsen  (i),  and  the  quadrupole  <9-  and  5-branch  ( 1-0)  band  measurements  of  Reuter  et  ai.  (6).  The  Q- 
and  5-branch  ( 1-0)  band  Raman  measurements  of  Bendtsen  (J)  were  excluded  from  the  fit  because  they 
show  a  bias  of  about  -f 0.009  cm~'  with  respect  to  our  measured  line  positions  and  our  fitting  results:  this 
offset  was  noted  previously  by  Reuter  et  al.  (6).  Similarly,  the  calibrated  coherent  anti-Stokes  Raman 
(CARS)  ( I-O)  band  ^-branch  measurements  of  Gilson  et  a!.  (5)  were  excluded  from  the  fit  because  of  a 
calibration  bias  in  the  data  (also  about  -f0.009  cm“’ ).  No  obvious  offset  was  noted  between  our  measurements 
and  the  0-branch  (1-0)  band  positions  from  Ref  (i)  so  we  retained  those  data  in  the  analysis. 

Table  I  presents  the  two  sets  of  Dunham  coefficients  obtained  from  the  least-squares  analysis.  Set  A  was 
obtained  by  fitting  the  coefficients  Tjo,  Xoi,  Xn,  and  X02  with  all  other  Dunham  coefficients  constrained  to 
zero.  For  set  B.  the  coefficient  X12  was  fitted  in  addition  to  the  same  four  Dunham  coefficients.  Also,  X2i 
was  fixed  at  the  value  -3.28  X  10~^  as  determined  by  Lofthus  and  Krupenie  (75),  and  X03  was  set  to  4.78 
X  10“’%  a  value  estimated  from  Dunham's  theory  (77)  by  substituting  the  Xqi  and  Xn  values  from  set  A 
and  the  Xjo  value  of  2358.5665  cm"’  from  Ref  (5)  into  Eq.  (2)  of  Reuter  et  al.  (6).  The  X,o  value  from 
our  analysis  was  not  used  in  these  calculations  because  only  two  vibrational  states  were  studied  in  this  work 


TABLE  II 

Experimental  Line  Positions^  for  ’"^NL 


V 

'  J' 

V"  J' 

.Measured 
'  Position 

R 

U 

V 

'  J' 

V 

Measured 

Position 

R 

U 

0 

2 

0  0 

11.9344^^ 

-2.9 

10.0 

1 

9 

0 

11 

2244.8345" 

-5.2 

3.0 

0 

3 

0  1 

19.8941*’ 

-0.9 

10.0 

1 

8 

0 

10 

2253.0967® 

-0.1 

3.0 

0 

4 

0  2 

27.8517*’ 

-0.3 

10.0 

1 

7 

0 

9 

2261.3203* 

-1.5 

10.0 

0 

5 

0  3 

35.8084*’ 

0.4 

10.0 

1 

6 

0 

8 

2269.5146® 

0.2 

3.0 

0 

6 

0  4 

43.7629*’ 

0.0 

10.0 

1 

5 

0 

7 

2277.6737* 

-0.6 

10.0 

0 

7 

0  5 

51.7163*" 

0.2 

10.0 

1 

4 

0 

6 

2285.8056* 

4.3 

10.0 

0 

8 

0  6 

59.6658*’ 

-1.8 

10.0 

1 

3 

0 

5 

2293.8986* 

3.5 

10.0 

0 

9 

0  7 

67.6162*" 

-0.9 

10.0 

1 

2 

0 

4 

2301.9561* 

0.7 

10.0 

0 

10 

0  8 

75.5643** 

0.1 

10.0 

1 

1 

0 

3 

2309.9851* 

3.2 

10.0 

0 

11 

0  9 

83.5061^ 

-2.5 

10.0 

1 

0 

0 

2 

2317.9750* 

0.7 

10.0 

0 

12 

0  10 

91.4513*" 

1.1 

10.0 

1 

6 

0 

4 

2372.9433 

-1.6 

2.4 

0 

13 

0  11 

99.3881*" 

-0.5 

10.0 

1 

8 

0 

6 

2388.3288 

0.2 

0.2 

0 

14 

0  12 

107.3242^ 

0.7 

10.0 

1 

9 

0 

7 

2395.9652 

-0.2 

0.2 

0 

15 

0  13 

115.2546^ 

-0.1 

10.0 

1 

10 

0 

8 

2403.5651 

0.0 

0.2 

0 

16 

0  14 

123.1839^ 

2.0 

10.0 

1 

11 

0 

9 

2411.1273 

-0.0 

0.3 

0 

17 

0  15 

131.1061^ 

1-2 

10.0 

1 

12 

0 

10 

2418.6519 

-0.1 

0.2 

0 

18 

0  16 

139.0205** 

-2.8 

10.0 

1 

13 

0 

11 

2426.1389 

0.2 

0.2 

0 

19 

0  17 

146.9378* 

1.0 

10.0 

1 

14 

0 

12 

2433.5871 

-0.1 

0.2 

0 

20 

0  18 

154.8418*" 

-3.5 

10.0 

1 

15 

0 

13 

2440.9973 

0.0 

0.2 

1 

17 

0  19 

2177.6601 

1.1 

1.0 

1 

16 

0 

14 

2448.3709® 

2.4 

3.0 

1 

16 

0  18 

2186.1621 

-1.2 

0.8 

1 

17 

0 

15 

2455.7008 

0.1 

0.2 

1 

15 

0  17 

2194.6277* 

-9.9 

30.0 

1 

18 

0 

16 

2462.9934 

-0.2 

0.3 

1 

14 

0  16 

2203.0818 

0.1 

0,4 

1 

19 

0 

17 

2470.2474® 

0.5 

7.0 

1 

13 

0  15 

2211.4956 

0.2 

0.5 

1 

20 

0 

18 

2477.4619® 

1.6 

3.0 

1 

12 

0  14 

2219.8781 

-0.2 

0.5 

1 

21 

0 

19 

2484.6349 

1.4 

1.9 

1 

11  • 

0  13 

2228.2309 

0.7 

0.9 

1 

22 

0 

20 

2491.7658 

-0.5 

1.3 

1 

10 

0  12 

2236. 5481" 

-2.6 

3.0 

_ 

_ 

a)  Measured  line  positions  are  in  R  is  the  residual  (measured  minus 

calculated)  in  10"^  cm'^  obtained  with  the  set  B  constants  given  in  Table  I, 
and  U  is  the  estimated  uncertainty  in  10“^  cm~^.  Uncertainties  for 

transitions  from  Refs.  (3,6)  are  taken  from  Table  I  of  Ref. (6). 
Uncertainties  from  the  present  work  have  been  estimated  as  described  in  the 
text . 


b) 

c) 


These  measurements 
These  measurements 


are  taken  from  Ref. 
are  taken  from  Ref. 


(3). 

(6). 
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and  all  higher-order  pure  vibrational  Dunham  coefficients  were  constrained  to  values  of  zero.  This  procedure 
produ  es  a  fitted  coefficient  that  is  an  underestimate  of  y,„.  In  Table  1.  we  also  give  the  weighted  fitting 

eJrr(A'Z  w.ar/rit  I  where  A' is  the  number  of  lines.  is  the  number  of  degrees  of  freedom,  and 

0  is  the  observed  minus  calculated  wavenumber  for  line  1(6).  a-  ■  .  ,  „mk, 

'  Examination  of  the  values  in  Table  I  indicates  that  the  uncertainties  in  the  case  A  coefficien  s  sh^tly 
smaller  than  those  of  case  B.  The  case  A  and  case  B  values  are  consistent  among  themselves,  but  the 
coefficients  from  the  two  solutions  are  sometimes  slightly  outside  the  range  of  the  calculated  uncertainties. 
Table  I  also  lists  previouslv  reported  values  for  the  Dunham  coefficients.  The  accuracies  of  >o,,  To:.  T.o, 
IS  y  ,  hSe  Sn  improv;d  by  about  a  factor  of  3  compared  to  the  results  in  Ref  (6).  No  eiror  estimates 
were  given  for  the  Dunham  coefficients  listed  in  Ref.  (5),  but  the 

with  those  results.  The  IV  coefficient  determined  in  case  B  has  an  uncertainty  of -50%.  Within  this  un- 
certaintv  the  measured  value  of -1.01  X  iQ-*  agrees  with  a  value  of -7.9  X  10  calculated  froin  EM 
of  Reuter  ei  al.  (6).  the  case  B  values  of  y„  and  To,,  and  the  values  of  y.o  and  y^o  reported  by  Gilson 

"'^  -fSle  11  presents  the  measurement  results  and  the  residuals  of  the  case  B  weighted  least-squares  best  fit 
to  the  dau  The  case  A  results,  which  are  not  reported,  are  very  similar.  For  most  lines,  the  residual  ( me^ured 
mils  cLlate^^^^^^^  is  less  than  the  estimated  uncertainty  .  The  residuals  of  both  the  far  infrared  and 
m  data  show  no  obvious  trend  with  J.  Reuter  «  al  (6)  noted  that  the  pure  rotational  line^  showed  a 
Slight  negative  bias  in  the  residuals  of  their  fit,  but  this  offset  does  not  appear  m  the  present  results  nor  is  i 
obvious  in  ground  stale  combination  differences  computed  from  the  stratospheric  infrared  line 

We  examined  published  0.0033*cm-’  resolution  stratosphenc  emission  spectra  for  the  ^ 

ground  state  f  =  0-0  quadrupole  lines  of  N,.  The  7-40  cm'^  interval  covered  m  the  atlas  J 

ll  {]9  m  and  the  few  small  higher  wavenumber  intemls  in  Refs.  (21-23)  show  only  the  S{  1 )  an  (  ) 
fin«  at  19  895  and  43.763  cm-  in  regions  nearly  free  of  interfering  tellunc  features  (see  Fig.  13  of  Ref 
(19)  and  Fig.  2  of  Ref  (23)).  The  absence  of  any  significant  features  at  the  predicted  positions  o 
transitions  indicates  that  the  N.  pure  rotational  quadrupole  band  is  extremely  weak. 
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Analysis  of  Atmospheric  Trace  Constituents 
from  High  Resolution  Infrared  Balloon-Borne  and 
Ground-Based  Solar  Absorption  Spectra 
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AB.STRACT 

Recent  results  and  ongoing  studies  of  high  resolution  solar  absorption  spectra  wll  be  presented.  The  analysis  of  these 
spectra  is  aimed  at  the  identification  and  quantification  of  trace  constituents  important  in  atmospheric  chemistry  of  the 
stratosphere  and  upper  troposphere. 

Analysis  of  balloon-borne  and  ground-based  spectra  obtained  at  0.0025  cm'*  covering  the  700-2200  cm'*  interval  wll  be 
presented.  Results  from  ground-based  0.02  cm*‘  solar  spectra,  from  several  locations  such  as  Denver,  South  Pole,  M.  Loa, 
and  New  Zealand  will  also  be  shown.  The  0.0025  cm'*  spectra  show  many  new  spectroscopic  features.  The  analysis  of 
these  spectra,  along  vdth  corresponding  laboratory  spectra,  improves  the  spectral  line  parameters,  and  thus  the  accuracy 
of  trace  constituents  quantification.  The  combination  of  the  recent  balloon  flights,  with  earlier  flights  data  since  1978  at 
0.02  cm'*  resolution,  provides  trends  analysis  of  several  stratospheric  trace  species.  Results  for  COFj.  F22,  SF4  and  other 
species  will  be  presented.  Analysis  of  several  ground-based  solar  spectra  provides  trends  for  HCl,  HF  and  other  species. 

The  retrieval  methods  used  for  total  column  density  and  altitude  distribution  for  both  ground-based  and  balloon-borne 
spectra  will  be  presented.  These  are  extended  for  the  analysis  of  the  ground-  based  spectra  to  be  obtained  by  the  high 
resolution  interferometers  of  the  Network  for  Detection  of  Stratospheric  Change  (NDSC).  Progress  of  the  University  of 
Denver  studies  for  the  NDSC  will  be  presented.  This  will  include  intercomparison  of  solar  spectra  and  trace  gases 
retrievals  obtained  from  simultaneous  scans  by  the  high  resolution  (0.0025  cm'*)  interferometers  of  BRUKER  and 
BOMEM. 


1  INTRODUCTION 

The  technique  of  infrared  solar  spectra  has  been  used  for  many  years  to  obtain  information  about  the  column 
abundance  and  vertical  distribution  of  atmospheric  gases.  To  maximize  retrieval  of  this  information,  nonlinear 
multiparameters  spectral  least  squares  fitting  methods  have  been  developed,  and  since  1980  became  an  important  tool 
in  our  research.*"*  Subsequent  studies  (see  below),  conducted  with  higher  resolution  and  signal  to  noise  spectra,  allowed 
to  further  enhance  these  analysis  methods. 

The  spectral  least  square  programs  can  fit  simultaneously  several  spectral  regions  of  several  spectral  scans,  with  three 
kinds  of  parameters:  (i)  instrumental  parameters,  such  as  resolution  and  background  contours  with  channel  and  phase 
corrections,  (ii)  atmospheric  (or  laboratory)  parameters  such  as  the  temperature  and  gas  amounts  m  isothermal  layers 
and  (iii)  spectral  line  parameters  such  as  line  positions,  intensities  and  half  widths.  The  programs  allow  the  initial 
parameters  to  remain  fixed  or  vary  within  a  constrained  range  of  a  priori  values.  The  fitted  parameters  include  the 
scaling  of  the  initial  mixing  ratio  distribution  functions  which  determine  the  gas  amount. 

The  atmospheric  spectroscopy  group  at  the  University  of  Denver  has  been  studying  atmospheric  infrared  absorption 
with  several  types  of  balloon-borne,  aircraft  and  ground-based  spectrometers  for  many  years.  The  instrumentation  used 
in  recent  experiments  includes  0.002-0.003  cm'*  (unapodized)  resolution  FTS  systems^  and  0.0005  cm'  resolution  tunable 
diode  LHS  system*  for  atmospheric  solar  absorption  spectra.  Ground-based  measurements  have  been  made  from  various 
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locations,  c.g.  Denver,  New  Mcmco,  j  locations  i!  the  U.S.  and  Europe.  The  field  measurements 

the  last  few  years.  Balloon  fl.ghts  were  th  the  field  spectrometers.  The  analysis  of  these  spectra 

arc  supplemented  by  extensive  laboratory  spectra,  ^  ,  atmospheric  trace  gases,  some  monitored  stnee  1980, 

has  been  very  successful  in  the  identification  an  qua  ^  Analysis  of  the  data  includes  detailed  atlases  with 

and  upgrading  the  theoretical  capabUiU^  of  atmospheric  spe^^  Secular  spectroscopic  constants  and  fine 

complete  identifications  of  the  individual  ^  ^  ^^^35  sections  simulations  of  emiss^^  and  absorption 

parameters  analysis  of  selected  absorption  parameters.  These  show  good 

spectra,  and  spectral  least-squares  fitting  using  the  most  c  spectral  features  in  the  stratospheric  spectrum, 

agreement  with  the  data  in  many  spectral  leading  to  advances  in  atmospheric  spectroscopy 

i*  ,hU  paper,  rumples  of  o«r  rocoo.  «  f™”  <>■“ 

wll  be  discussed.  Quantification  of  several  trace  gases  since  l^BU  win  o 

methods. 

RESULTS 


7 1  Ralloon-hnrne  spectra 

We  have  slodicd  Ihe  new  infrared  spectroseopic  strafosphcric  species 

„i,h  a  0.002  em"  resoWon  in.erferon.eter  system  "»»  HO.NO.  NO.,  COF.,  SF,  and 

arc  measurable  m  these  spectra.  The  mam  mo  spectra  and  laboratory  spectra  reveal  new  details  of  the 

CHC«F,.  The  stratospheric  solar  spectra,  ground-based  so  P  >  ^  .o^j^cmistry  and  the  heat  balance  of  the 

absorption  by  these  molecules  which  JJevbuW  unidentified  features  to  these  gases  in  the  stratospheric 

Many  of  these  features  were  not  accounted  for  in  SltsS  irdrall.^^^ 

dSlneThc  isotopic  inlios  of -‘O^O'-O  and -W-O^O  in 

The  cxcciicnt  agreement  between  the  h«hiy  prcc.se  Thc^^^^  resolution  of  the 

evidence  for  absorption  by  the  HO, NO.  8027  cm J,  ^ 

flights  measurements  reveals  additional  detai^  o  band,  and 

including  the  identification  of  overlapping  HNO3  ™  both  laboratory  and  stratospheric  spectra  New 

recent  analysis  of  the  HN03.g  and.  ,3,,  been  identified  in  the 

l^ictf^cr^^^^^^  an^-study  ^.h  new  iine  parameters  and  cross.ections. 

The  5-7  Jim  spectra  are  also  being  studied  for  in  the  6  urn  region  show 

bands,  and  H^O  .3  bands  are  now  used  for  new  etrieyals  of  HNO,^  The  ^  z  ^  ^  3ii„,,.3 

"r  iSe  i-e'Xr 

The  long  time  span  covered  by  medium  and  high  resolution  recent  0.002  enf*'  spectra 

considerable  interest  in  studying  trends  of  atmospheric  trace  g  ^  ^bc  last  decade.  Exponential 


22  nrotind-hn.scd  sncctra 


We  have  performed  quantitative 


itatlve  analysis  of  the  absorption  by  a  number  of  minor  and  trace  constituents  in  0.02  cm" 
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rcsoludon  solar  absorplio.  spectra  geographical  locatioas,-  In  ad«on  to  td 

include  the  total  vertical  column  amounts  of  HjO  (and  isotopes),  COj,  O3,  NjO,  4,  *  k  ^  losn 

HNO.  from  M.  Lea,  Hawaii,  in  February  1987-,  O.  N.O.  HNO.,  CXl.  CH  and  CF,C..  m  Dec^ber  mO  d 
™cmto-Dcccmbcr  1986  from  the  South  Pole;  HCI,  NO,  NO,  and  QH.  m  December  1986  from  S  ^ 
HO  HNO  O  CIONO,  and  NO,  from  McMurdo,  Antarctica  during  the  spring  1987  ozone  hole.  I”  the  M,  Loa 
mSu™ ^n'ts  U  ^  found  that  the  average  tropospheric  concentrations  deduced  lor  CO,  N,0  and  CH,  are  very 
consistent  with  correlative  NOAA  GMCC  surface  data. 

Of  oarticular  interest  U  the  observed  increase  in  the  CF.Ct,  column  amount  above  the  S.  Pole  >l>0  «  year  period, 
co°LCd“g  to  (3.6  ±  Zl)%  year-.  The  column  amount  of  HNO,  for  both  1980  and  1986  were  not  sigmlicantly 
aS  oler  data  sets,  sich  L  from  the  National  Solar  Observatory  (NSO)  at  Knt  Peak  were  used  for  mom  oring 
mends  of  HF  and  HC(  fro^n  0.02AI.005  cm-  spectra  obtained  during  the  1977-1990  period.  Unear  increase  rates  of  (10.7 
±  1.1)%  year*  and  (4.9  ±  0.7)%yeaf‘  have  been  obtained  for  HF  and  HC<  respectively. 

Wc  have  also  been  conducting  0.002-0.003  cm''  resolution  solar  absorption  ground-based  measurements  for  periodic 
ob^^lions  and  monitoring  of  the  atmosphere  from  Denver.  The  initial  measurements  were  made  w^J  the  Mloon- 
bornc  BOMEM  FTS  system.^  Typical  results  arc  shown  in  Figures  2-5.  More  rcccnll)^  dcchw  c 
BRUKER  FTS  have  been  used,  one  of  which  will  be  moved  to  Mauna  Loa  as  soon  as  the  NDSC  station  becomes 

operational.  Preliminary  results  are  shown  in  Figures  6-10. 

The  analysis  of  the  ground-based  spectra  shows  that  the  trace  gases  of  interest  to  the  NDSC  can  be 
from  selected  spectral  intervals  of  multiscans  measurements,  with  temperatures  and  geometry  an  i 
Xmance  moLted  from  reference  gaaca.  It  Is  apparent  that  no  more  than  scale  height  vertieal  ^ 

be  expected  for  the  altitude  distribution  of  the  trace  gases,  and  in  some  cases  only  the  tropospheric  and  f 

components  could  be  separated.  The  translation  of  the  statistical  errors  of  the  spectral  least  squares  fitting  to  abso 
errors  estimates  from  a  priori  constraints  and  accuracies  requires  further  studies. 

H  rONCLUSIONS 

Hich  resolution  stratospheric  spectra  show  many  new  spectral  features  not  observed  in  prewous  lower  resolution 
spectfa.  This  requires  new  analysis  of  laboratory  and  atmospheric  spectra  and  improved  theoretical 
molecular  bands  of  atmospheric  interest.  Atmospheric  spectra  remain  a  powerful  source  for  observing  and  monitoring 

our  changing  atmosphere. 

Retrieval  methods  based  on  spectral  least  squares  arc  best  suited  for  spectroscopic  quantification,  but  further 
development  are  desired. 
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rig.l^amplc  spectra  from  the  new  atlas  of  high  resolution  (0.002  cm'’)  stratospheric  IR  absorption  spectra  in  the  1250- 
1252  cm’*.  The  spectral  region  contains  severjJ  manifolds  of  the  COFj  1/4  band. 


OZONE  DU  JAN  12  1909  H=1.6KM 


Hg.2JJonlincar  tet-squares  mulli-sam  fit  of  O,  with  improved  line  paramclra  in  the  W 

obtained  from  Denver,  Colorado,  daring  saimet  on  Janaar,  12,  19^,  wrth  a  nrodrM  BOMEM  model  DA3.002 
interferometer  used  for  DU  balloon  flights.  The  total  column  of  the  fitted  profile  is  7.46x10  raolec/cm  . 
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OZONE  DU  JAN  12  1989  H=1.6KM 


ng.3.Nonlincar  least-squares  multi-scan  fit  of  O,  with  improved  line  parameters  m  the  1102  cm  region,  to  solar  spe^a 
obtained  from  Denver,  Colorado,  during  sunset  on  January  12,  1989,  with  a  modified  BOMEM  model  DA3.002 
interferometer  used  for  DU  balloon  flights.  The  total  column  of  the  fitted  profile  is  733x10  molec/cm  . 


C02  DU  JAN  12  1909  M^I.SKM 


rig.4.Nonlincar  least-squares  single  scan  fit  of  COj  in  the  936  cm*'  region  from  solar  spectra  obtained  from  Denver 
Colorado  during  sunset  on  January  12,  1989,  with  a  modified  BOMEM  model  DA3.002  interferometer  used  for  DU 
balloon  flights. 
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HN03  DU  JAN  12  1989  II=1.6KM 


Fig.S.Nonlincar  least-squares  multi-scan  fit  of  HNO3  with  improved  line  parameters  in  the  868  cm  region,  to  solar  spectra 
obtained  from  Denver,  Colorado,  during  sunset  on  January  12,  1989,  with  a  modified  BOMEM  model  DA3.002 
interferometer  used  for  DU  balloon  flights.  The  total  column  of  the  fitted  profile  is  1.35x10*'^  molcc/cml 


COADD  SOL  BRUK  COADD  8  SOL  BOU 


Fig.6.Solar  spectra  obtained  from  Denver,  Colorado,  during  preliminary  testing  of  the  two  interferometer  systems 
BOMEM  DA8.002  and  BROKER  IFS  120  HR  in  the  COj  935-940  cm'^  region. 
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WAVENUMBF.n  (cm'*  ) 

Fig.7.SoIar  spectra  obtained  from  Denver,  Colorado,  during  preliminary  testing  of  the  two  interferometer  systems 
BOMEM  DA8.002  and  BRUKER  IPS  120  HR  in  the  HNO3  865-870  cm'^  region. 
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Fig.8.Solar  spectra  obtained  from  Denver,  Colorado,  during  preliminary  testing  of  the  two  interferometer  systems 
BOMEM  DA8.002  and  BRUKER  IPS  120  HR  in  the  O3  996-998  cm''  region. 
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APPENDIX  9 

a  reprint  from  Applied  Optics 


Improved  line  parameters  for  ozone  bands  in 
the  10-)Lim  spectral  region 


Jean-Marie  Flaud,  Claude  Camy-Peyret,  Curtis  P.  Rinsland,  V.  Malathy  Devi,  Mary  Ann  H.  Smith,  and 
Aaron  Goldman 


A  complete  update  of  spectroscopic  line  parameters  for  the  bands  of  ozone  is  reported.  The  listing 

contains  calculated  positions,  intensities,  lower  state  energies,  and  air-  and  self-broadened  halfwidths  of  more 
than  53,000  lines.  The  results  have  been  generated  using  improved  spectroscopic  parameters  obtained  in  a 
number  of  recent  high  resolution  laboratory  studies.  A  total  of  eighteen  bands  of  ^^63  (sixteen  hot  bands  plus 
the  vi  and  V3  fundamentals)  are  included  along  with  the  and  V3  fundamentals  of  both  and 

leQiso^^O.  As  shown  by  comparisons  of  line-by-line  simulations  with  0.003-cm"^  resolution  balloon-borne 
stratospheric  solar  spectra,  the  new  parameters  greatly  improve  the  accuracy  of  atmospheric  calculations  in 
the  10’fj.m  region,  especially  for  the  isotopic  and  lines. 


The  need  for  improvements  in  the  accuracy  of  re¬ 
mote  sounding  measurements  of  ozone  in  the  earth’s 
atmosphere  has  motivated  a  number  of  recent  studies 
devoted  to  the  measurement  and  calculation  of  the 
spectrum  of  this  molecule.  In  the  IR,  the  strongest 
bands  are  at  10  fim,  where  absorption  is  produced 
primarily  by  the  and  vi  fundamentals  of  Hot 

bands  of  and  the  ps  and  Pi  bands  of  and 

160180^^0  also  contribute  significantly  to  the  total  ab¬ 
sorption  by  ozone.^  The  IO-mdi  region  has  been  used 
in  a  number  of  atmospheric  studies  of  ozone,  for  exam¬ 
ple,  the  measurement  of  stratospheric  O3  profiles  by 
the  Limb  Infrared  Monitor  of  the  Stratosphere 
(LIMS)  instrument  on  the  Nimbus  7  satellite, mea¬ 
surements  of  ozone  total  columns  inside  the  Antarctic 
ozone  hole  using  high  resolution  ground  based  solar 
spectroscopy,^  and  the  measurement  of  ozone  isotope 
ratios  from  ground  based  and  balloon-borne  high  reso¬ 
lution  solar  spectra. Accurate  knowledge  of  the  ab¬ 
sorption  properties  of  the  10-fim  ozone  bands  is  also 
needed  for  climate  models  to  estimate  future  changes 
in  atmospheric  and  surface  temperatures  produced  by 
the  greenhouse  effect, 

In  the  present  paper,  we  report  the  calculation  of 
improved  line  parameters  for  the  ozone  bands  in  the 
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10-)um  region.  In  addition  to  revising  the  eight  bands 
on  the  1986  HITRAN^^  and  1984  GEISA^^  line  param¬ 
eter  compilations,  the  calculations  have  been  extended 
to  include  fourteen  additional  bands.  For  example, 
only  four  hot  bands  of  the  Pi  +  P2'-  p2,  1^2  +  —  p2, 

I'l  +  I's  -  z'l,  and  2z^3  -  were  included  previously,^^-^^ 
with  upper  state  energy  levels  derived  from  spectra 
recorded  at  0.02-cm”^  resolution. The  new  pa¬ 
rameters  listing  includes  a  total  of  sixteen  hot  bands  of 
^^03  based  on  extensive  analysis  of  positions  and  inten¬ 
sities  from  0.005-  to  0.01-cm“^  resolution  laboratory 
spectra.  Also,  for  the  first  time,  accurate  parameters 
for  the  and  vz  bands  of  and  are 

included.  The  pz  bands  of  both  isotopic  species  ap¬ 
peared  in  the  1986  HITRAN  and  1984  GEISA  compi¬ 
lations, but  the  crude  approximations  used  to  gen¬ 
erate  the  line  list  make  the  data  unreliable  for  high 
resolution  studies.  All  the  new  results  have  been  de¬ 
rived  from  laboratory  absorption  spectra  of  ozone  re¬ 
corded  with  the  McMath  Fourier  transform  spectrom¬ 
eter  operated  on  Kitt  Peak  by  the  National  Solar 
Observatory.^^““^ 

A  summary  of  the  new  parameters  is  given  in  Table  L 
The  integrated  band  intensities  Sband  have  been  evalu¬ 
ated  without  a  minimum  intensity  cutoff  applied  in  the 
computer  code.  This  procedure  gives  the  most  accu¬ 
rate  estimate  for  the  integrated  band  intensity  but 
includes  many  very  weak  lines  not  readily  observable 
in  long  path  atmospheric  or  laboratory  spectra.  For 
this  reason,  the  line  listing  was  generated  with  mini¬ 
mum  line  intensity  cutoffs  of  1.0  X  lO^^s  cm“Vmole" 
cule  cm”^  at  296  K  for  the  bands  and  5.0  X  10“^^ 
cm’Vmolecule  cm^^  at  296  K  for  the  bands  of  the  two 
^®0-monosubstituted  ozone  isotopes.  The  sum  of  in¬ 
tensities  in  this  list,  denoted  Snsu  is  also  reported  in 
Table  I.  The  calculations  include  a  scaling  of  the 
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Table  1. 

SufTimary  of  O3  Bands  Included  in  the  Calculations 

Vibrational 

Calculated 

Nmaber 

C  ►  ^ 

Band  Center 

Intenst. 

ties- 

Upper 

Lover 

Isotope* 

(cm"^) 

Sa.t'’ 

Lines 

012 

110 

666 

929.8447 

4.16x10"“ 

3.29x10'“ 

164 

002 

100 

666 

954 .7535 

3.33x10'“ 

3.27x10'^* 

951 

lU 

110 

666 

988.9772 

2.12x10'“ 

2,10x10'=* 

1342 

012 

oil 

666 

999.5841 

6.01x10"^^ 

6.00x10'“ 

1603 

101 

100 

666 

1007.6470 

6.25x10'“ 

6.25x10'“ 

2546 

001 

000 

686 

1008.4528 

-2.65x10'“ 

2.65x10'“ 

2095 

021 

020 

666 

1008.6618 

0 

0 

1.50x10*“ 

I5U 

002 

001 

666 

1015.8068 

1.68x10'“ 

1.68x10'“ 

3137 

Oil 

010 

666 

1025.5914 

4.64x10"“ 

4.64x10'“ 

3883 

001 

000 

668 

1028.1120 

5.12x10"” 

5.12x10'“ 

3827 

001 

000 

666 

1042.0840 

1.41x10'“ 

1.41x10*“ 

7224 

111 

oil 

666 

1058.7166 

1.22x10'“ 

6.38x10"“ 

252 

101 

001 

666 

1068.7003 

5.31x10"“^ 

5.24x10*“ 

2558 

100 

000 

686 

1074.3076 

4.54x10'“ 

4.54x10'“ 

838 

120 

020 

666 

1087.3041 

2.72x10'“ 

2.21x10'“ 

820 

210 

110 

666 

1089.9162 

8.85x10'“ 

3.33x10'“ 

177 

100 

000 

668 

1090.3541 

4.07x10'^^ 

4.07x10'=^ 

4425 

ILO 

010 

666 

1095,3308 

1.27x10'“ 

1.27x10*“ 

3885 

200 

100 

666 

1098.0179 

3.37x10'“ 

3.30x10*=* 

2931 

100 

000 

666 

1103.1373 

5.40x10'“ 

5.40x10'“ 

6766 

200 

001 

666 

1159.0712 

8.42x10'“ 

8.34x10'=^ 

2057 

210 

oil 

666 

1159.6556 
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intensities  vithout  a  minimum  intensity  cutoff. 

Sum  of  the  intensities  of  all  lines  in  linelisC,  The  linelist  was 
calculated  with  minimum  intensity  cutoffs.  See  text. 


intensitiss  by  a  multiplicativa  factor  to  convert  the 
intensities  from  values  for  a  pure  isotopic  species  to 
values  for  a  sample  in  natural  isotopic  molecular  abun¬ 
dance.  For  this  purpose,  we  adopted  the  isotopic  ra¬ 
tios  on  the  HITRAN  compilationio  (isOg  =  0.9928, 
160160180  =  0.0040,  and  i^oisoi^O  =  0.0020).  The 
differences  between  Sband  “^ust  significant  only 
for  the  very  weak  bands. 

It  should  be  noted  that  recent  studies  (see  Ref.  6  ^d 
references  cited  therein)  have  indicated  that  the  iso¬ 
topic  abundances  of  the  two  i®0-mono8ubstituted 
ozone  isotopes  are  enhanced  in  the  stratosphere  com¬ 
pared  to  the  expected  normal  isotopic  values.  There¬ 
fore,  atmospheric  simulations  generated  with  the  new 
line  parameters,  which  incorporate  the  natoal  abun¬ 
dance  factors  used  in  the  HITRAN  compilation,^®  may 
underestimate  the  absorption  produced  by  the 
160160^^0  and  lines.  The  enrichment  effect 

can  be  simulated  by  scaling  the  and 

160180160  line  intensities  by  the  appropriate  isotopic 
enhancement  factors  in  each  atmospheric  layer. 

The  calculations  of  the  positions  and  intensities  of 
the  vz  and  Pi  fundamentals  of  and  the  two  i^O- 
monosubstituted  ozone  isotopes  have  been  described 
previously. Since  only  relative  intensities  could 
be  determined  from  the  laboratory  spectra,  it  is  impor¬ 
tant  to  recall  the  assumptions  used  in  scaling  the  rela¬ 
tive  intensities  to  obtain  absolute  values.  For  ^®03,  a 
value  of  ”  "0.2662  D  was  assumed,  and  for 

160180160  the  intensities  were  computed  by  transfer¬ 
ring  the  dipole  moment  of  to  i60^®0^®0.  The 
160160180  intensities  have  been  calculated  from  rela¬ 


tive  intensities  of  i60i®0i®0  and  isqisoisO  lines  mea¬ 
sured  in  the  experimental  spectra.^®  The  amount  of 
160180160  in  the  cell  was  first  determined  by  compar¬ 
ing  the  measured  relative  intensities  of  the  i6q‘80i60 
lines  with  the  computed  absolute  values.^"  This 
amount  was  then  multiplied  by  a  factor  of  2  to  account 
for  the  difference  in  the  relative  concentrations  of 
160160180  and  i60i®0i®0. 

The  line  positions  for  the  hot  bands  V2  +  vz~  *'2  and 
vi  +  VO  ~  "2  of  ^^Og  were  computed  from  improved 
energy  levels  for  the  upper  and  lower  states  of  both 
bands.  The  (Oil)  and  (110)  energy  levels  were  derived 
from  analysis  of  the  vi  +  vg  and  pz  "I"  *'3  bands  in  the  5.7- 
/im  region.^®  The  energy  levels  for  the  (010)  state  were 
calculated^®  from  a  refit  to  line  positions  measured  in 
the  P2  and  2^2  -  p2  bands.  The  intensities  of  these 
hot  bands  were  computed  assuming  transition  mo¬ 
ments  derived  from  those  of  the  I'l  and  pz  bands. 

The  line  positions  for  the  {(002), (101), (200))- 
{(100),(001)}  hot  band  system  of  ^^Og  were  computed 
from  improved  energy  levels  derived  for  the  2j'3,  I'l  -f-  vg, 
and  2pi  bands  in  the  4.8-Mni  region-^  and  the  funda¬ 
mental  bands  at  10  pm.^®  As  mentioned  previously,^! 
the  observed  energies  of  the  (101)  [19  2  17]  and  (lOH 
[45  3  42]  rovibrational  levels  are  displaced  from  their 
calculated  values  probably  because  of  a  resonance  with 
an  unidentified  nearby  level.  The  calculated  energies 
of  both  levels  were  replaced  with  their  observed  values 
in  generating  the  line  list.  The  intensities  were  com¬ 
puted  from  transition  moments  derived  from  those  for 
the  Pz  and  Pi  bands  and  assuming  the  higher-order 
dipole  moment  derivatives  are  equal  to  zero  [except  for 
the  transition  moment  of  (002)-(001)  which  was  deter¬ 
mined  to  be  0.5141  X  lO"^  D  from  a  least-squares  fit  of 
measured  intensities]. 

Positions  for  the  lines  from  the  {(021), (120)}— (020) 
system  have  been  calculated  from  energy  levels  of  the 
2^2  +  vz  and  2^2  +  pi  bands  obtained  from  the  recent 
first  high  resolution  analysis  of  these  bands^^  and  the 
constants  for  the  (020)  state  reported  by  Flaud  etal.^ 
The  intensities  were  computed  assuming  transition 
moments  derived  from  those  of  the  vi  and  pz  bands. 

Although  the  transitions  are  rather  weak  at  atmo¬ 
spheric  temperatures  because  of  the  Boltzmann  factor, 
the  {(012),(111),(210)}-((110),(011)}  system  has  been 
included  to  complete  the  description  of  the  lO-^m 
region.  The  positions  were  calculated  from  upper 
state  energy  levels  determined  from  analysis  of  the  pz  + 
2:'3,  Pi  +  P2  +  vs,  and  2pi  +  pz  bands  in  the  3.6-Mm 
region^®  and  lower  state  energy  levels  determined  from 
analysis  of  the  Pi  +  pz  and  pz  4"  vz  bands  in  the  S-l-fim 
region.!®  The  intensity  calculation  has  been  per¬ 
formed  assuming  that  the  transition  moments  of 
(210)-(011)  and  (012)-(110)  are  equal  to  zero,  and  the 
other  ones  were  derived  from  those  of  the  vi  and  vz 
bands. 

The  accuracies  of  the  calculated  line  positions  and 
intensities  of  the  pz  and  Pi  bands  of  !®03,  !60!®0!®0,  and 
i60180!60  were  given  in  the  previous  papers.!®"!’ 
Also,  the  accuracies  of  the  calculated  positions  of  the 
!803  hot  band  lines  can  be  estimated  by  combining  the 
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uncertainties  of  the  calculated  upper  and  lower  energy 
levels  as  reported  previously,  A  precise  error 

budget  is  difficult  to  estimate  for  the  hot  band  line 
intensities  since  they  have  been  calculated  from  transi¬ 
tion  moment  expansions  derived  from  those  of  the 
corresponding  cold  bands;  generally,  we  expect  the 
intensity  errors  of  ^^10%  with  the  accuracy  degrading 
with  increasing  quantum  numbers. 

As  an  interim  calculation,  we  have  included  air- 
broadened  and  self-broadened  halfwidths  based  on 
empirical  fits  to  recent  laboratory  measurements  of 


Fig.  1.  Comparison  of  a  0.003-cm“^  resolution  solar  occultation 
spectrum  and  a  simulation  of  the  data  generated  including  the  1986 
HITRAN  ozone  line  parameters'^  (upper  panel)  and  a  simulation  of 
the  data  generated  with  new  ozone  line  parameters  (lower  panel). 
In  each  panel,  the  upper  curve  is  the  line-by-line  simulation  which 
has  been  displaced  vertically  by  0.2  for  clarity.  The  lower  curve 
shows  a  solar  spectrum  obtained  at  sunset  6  June  1988  during  a 
balloon  flight  from  Palestine,  TX  conducted  by  the  University  of 
Denver  atmospheric  spectroscopy  group.  The  spectrum  (normal¬ 
ized  to  the  highest  intensity  in  the  interval)  was  recorded  with  a 
Michelson  interferometer  from  a  float  altitude  of  36.85  km  and  at  a 
solar  astronomical  zenith  angle  of  92.61®.  The  corresponding  re¬ 
fracted  tangent  height  is  30.27  km.  Arrows  and  capital  letters  in  the 
lower  panel  indicate  examples  of  ozone  spectral  features  produced 
by  different  vibration-rotation  bands.  The  ozone  isotopic  species 
and  vibrational  bands  corresponding  to  the  marked  features  are  A  = 
(002)-(100),B  =  leoiSQieo  (OOl)-(OOO),  C  =  (002)-(001),D 

=  160160180  (OOl)-(OOO),  E  =  ^603  (Oll)-(OIO),  F  =  (lOl)-(lOO), 
G  =  i6O3(100)-(000),H=  '^03  (001)-(000),and7=  (021)-(020). 

The  R28  line  of  the  00011-10001  band  of  at  980.9132  cm"^ 

(marked  by  a  solid  circle  in  both  panels)  is  the  only  significant  non- 
O3  absorption  feature  in  this  spectral  interval. 


ozone  halfwidths.-^-^^  The  following  expression  has 
been  adopted: 

^  <3o  +  ai</"  +  aoJ"^  +  i 

where  aL  is  the  Lorentz  broadening  coefficient  at  296  K 
(in  cm*"^  atm“^  units)  and  is  the  lower  state  total 
angular  momentum  quantum  number.  For  the  air- 
broadened  halfwidths,  the  coefficients  Co  =  0,08735,  Ui 
=  —0.0007447, 02  =  6.08  X  10“^,  and  03  =  0.0  have  been 
used.  On  average,  the  measured  air-broadened  half¬ 
widths  in  the  vj  are  ^^10%  larger  than  corre¬ 

sponding  values  on  the  1986  HITRAN  compila¬ 
tion. The  self-broadened  halfwidths  were 
calculated  with  the  coefficients  oq  =  0.11083,  oi  = 
-0.0016721,  02  =  5.0084  X  10*“^  and  03  =  -4.6741  X 
10“'^  for  J"  <  25.  For  J'"  >  25,  aL  =  0.093  has  been 
adopted.  Self-broadened  halfwidths  were  not  includ¬ 
ed  on  the  1986  HITRAN^^  and  1984  GEISA  compila¬ 
tions.^^ 

The  new  ozone  line  parameters  yield  considerable 
improvement  in  the  accuracy  of  high  resolution  atmo¬ 
spheric  transmission  and  emission  calculations  in  the 
lO-Atm  region,  especially  for  the  isotopic  and 

leoiSQieo  bands  and  the  hot  bands  of  ^^Os,  To  illus¬ 
trate  the  quality  of  the  new  line  parameters,  Figs.  1-3 
show  comparisons  of  synthetic  spectra  with  0.003- 
cm”^  resolution  stratospheric  solar  occultation  spectra 
recorded  with  a  balloon-borne  Fourier  transform  spec- 
trometer.^’2"’-®  In  all  cases,  the  agreement  between 
measurement  and  calculation  is  very  good,  even  for  the 
weakest  lines.  Figure  1  includes  a  comparison  of  the 
stratospheric  spectra  with  a  simulation  of  generated 
data  with  the  1986  HITRAN  compilation  ozone  line 
parameters.  The  limitations  of  the  previous  data  for 
high  resolution  spectroscopic  studies  are  obvious. 


Wavenumber  (cm 


Fig.  2.  Comparison  of  a  0.(X)3-cm“^  resolution  stratospheric  solar 
occultation  spectrum  (same  scan  as  shown  in  Fig.  1)  and  line-by-line 
simulations  with  the  new  ozone  line  parameters.  The  upper  spec¬ 
trum,  which  has  been  displaced  vertically  by  0.1  for  clarity,  shows 
the  absorption  computed  for  the  isotopic  species  only. 

The  feature  near  1090.36  cm”"^  is  the  Q-branch  of  the  i^i  band.  More 
than  300  lines  are  included  in  the  calculations.  The  two  lower 
spectra  show  a  comparison  of  the  measured  solar  spectrum  (solid 
line)  with  a  simulation  of  the  absorption  by  all  gases  (open  dia¬ 
monds).  The  absorption  in  this  region  is  due  almost  entirely  to  ^^63 
lines  and  the  lines  of  the  i^i-band  Q-branch. 
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Fig.  3.  Comparison  of  a  0.003-cm“‘  resolution  stratospheric  solar 
occultation  spectrum  (same  scan  as  shown  in  Figs.  1  and  2)  and  a 
line-by-line  simulation  including  the  new  ozone  line  parameters.  A 
number  of  features  of  the  (200)-(001)  band  of  *®03  are  identified. 


The  total  number  of  ozone  lines  in  the  new  10-;im 
line  list  is  53,735  compared  with  20,969  in  the  1986 
HITRAN  compilation.  10  The  substantial  increase  in 
the  number  of  lines  is  likely  to  affect  atmospheric 
retrievals  over  long  atmospheric  paths,  particularly 
limb-viewing  paths  through  the  lower  stratosphere 
where  the  ozone  abundance  is  a  maximum.  The  new 
line  list  also  includes  some  significant  changes  in  the 
distribution  of  the  intensities  within  the  stronger 
bands.  An  important  example  is  the  j^i  band  of  i^Os. 
Although  the  integrated  intensity  of  the  band  in  the 
new  line  list  is  only  slightly  different  (5.40  X  10  i® 
compared  with  5.56  X  10”i®  cm~Vmolecule  cm  ^  at  296 
K  in  the  1986  HITRAN  listi®),  the  intensities  of  indi¬ 
vidual  lines  selected  for  ozone  retrievals  from  lugh 
resolution  atmospheric  spectra  are  in  most  cases  high¬ 
er  than  on  the  1986  HITRAN  list,io  resulting  in  corre¬ 
spondingly  lower  inferred  ozone  amounts.^®’®®  The 
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Fig.  4.  Absorption  line  intensity  (cm-'/molecule  cm'^  at  296  K  units)  on  a  base  10  logarithm  vs  wavenumber  scale  for  each  ozone  band.  A 
plot  with  all  the  lines  is  presented  in  the  same  format  at  the  bottom  of  the  right  column. 
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retrievals  with  the  new  line  parameters  are  in  better 
agreement  with  correlative  ozone  measurements  than 
those  derived  with  the  earlier  line  list.-^’^^  Note  that 
changes  in  the  line  intensities  have  an  important 
effect  on  the  derivation  of  and 

160180160/16O3  isotopic  enhancement  ratios  from  high 
resolution  IR  spectra.  The  recent  analysis  of  0.003- 
cm^^  resolution  stratospheric  solar  absorption  spectra 
by  Goldman  et  aL^  determined  these  ratios  using  the 
line  parameters  reported  in  this  paper. 

To  provide  an  overview  of  the  structures  and  the 
relative  intensities  of  the  bands  in  the  new  line  list,  Fig. 
4  presents  plots  of  absorption  line  intensity  (cm^'V 
molecule  cm"-  at  296  K  units)  on  a  base  10  logarithm  vs 
wavenumber  scale  for  each  ozone  band.  The  bands 
are  ordered  by  increasing  band  center,  starting  at  the 
top  of  the  left  column  and  ending  at  the  bottom  of  the 
right  column.  The  final  plot  (lower  right)  contains  all 
the  lines  displayed  in  the  same  format.  The  irregular 
distribution  of  intensity  vs  wavenumber  within  the 
bands  results  from  the  numerous  resonances,  both 
smooth  and  localized,  occurring  between  the  interact¬ 
ing  states. 
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the  University  of  Denver  is  supported  under  coopera¬ 
tive  agreements  with  the  National  Aeronautics  and 
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References 

1.  A.  Goldman,  T.  G.  Kyle,  D.  G.  Murcray,  F.  H,  Murcray,  and  W.  J. 
Williams,  “Long  Path  Atmospheric  Ozone  Absorption  in  the  9- 
lO-Min  Region  Observed  from  a  Balloon-Borne  Spectrometer, 
Appl.  Opt,  9,  565-580  (1970). 

2.  E.  E.  Remsberg  et  aL^  “The  Validation  of  NIMBUS  7  LIMS 
Measurements  of  Ozone,”  J.  Geophys,  Res.  89,  5161-5178 
(1984). 

3.  S.  R.  Drayson  et  a/.,  “Spectroscopy  and  Transmittances  for  the 
LIMS  Experiment,”  J.  Geophys.  Res.  89, 5141-5146  (1984). 

4.  C.  B.  Farmer,  G.  C.  Toon,  P.  W.  Schaper,  J.-F.  Blavier,  and  L.  L. 
Lowes,  “Stratospheric  Trace  Gases  in  the  Spring  1986  Antarctic 
Atmosphere,”  Nature  London  329, 126-130  (1987). 

5.  C.  P.  Rinsland,  V.  Malathy  Devi,  J.-M.  Flaud,  C.  Camy-Peyret, 
M.  A.  H.  Smith,  and  G.  M.  Stokes,  “Identification  of  ^®0-Iso- 
topic  Lines  of  Ozone  in  Infrared  Ground-Based  Solar  Absorp¬ 
tion  Spectra,”  J.  Geophys.  Res.  90,  10,719-10,725  (1985). 

6.  A.  Goldman  et  a/.,  “Isotopic  Abundances  of  Stratospheric  Ozone 
from  Balloon-Borne  High  Resolution  Infrared  Solar  Spectra,”  J. 
Geophys.  Res,  94,  8467-8473  (1989). 

7.  V,  Ramanathan,  R,  J.  Cicerone,  H.  B.  Singh,  and  J.  T.  Kiehl, 
“Trace  Gas  Trends  and  Their  Potential  Role  in  Climate  Chan¬ 
ge,”  J.  Geophys.  Res.  90,  5547-5566  (1985). 

8.  R.  E.  Dickinson  and  R,  J,  Cicerone,  “Future  Global  Warming 
from  Atmospheric  Trace  Gases,”  Nature  London  319,  109-115 
(1986). 

9.  D.  P.  Kratz  and  R.  D.  Cess,  “Infrared  Radiation  Models  for 
Atmospheric  Ozone,”  J.  Geophys.  Res.  93,  7047-7054  (1988). 

10.  L.  S.  Rothman  et  ai,  “The  HITRAN  Database:  1986  Edition,” 
Appl.  Opt  26,  4058-4097  (1987). 

11.  N.  Husson  et  ai,  “The  GEISA  Spectroscopic  Line  Parameters 
Data  Bank  in  1984,”  Ann.  Geophys.  4, 185-190  (1986). 


12.  J.-M.  Flaud,  C.  Camy-Peyret,  and  L.  S.  Rothman,  “Improved 
Ozone  Line  Parameters  in  the  10-  and  4.8-Mni  Regions,”  Appl. 
Opt  19, 655(1980), 

13.  A.  Barbe,  C.  Secroun,  P.  Jouve,  C.  Camy-Peyret,  and  J.-M. 
Flaud,  “High  Resolution  Infrared  Spectrum  of  the  vo  +  and 

+  Vo  Bands  of  Ozone,”  J.  Mol.  Spectrosc.  75, 103-110  (1979). 

14.  J.-M.  Flaud,  C.  Camy-Peyret,  A.  Barbe,  C.  Secroun,  and  P. 
Jouve,  “Line  Positions  and  Intensities  for  the  21^3,  vi  +  vsy  and  2vi 
Bands  of  Ozone,”  J.  Mol.  Spectrosc.  80, 185-199  (1988). 

15.  J.-M.  Flaud,  C.  Camy-Peyret,  V.  Malathy  Devi,  C.  P,  Rinsland, 
and  M.  A.  H.  Smith,  “The  vi  and  vs  Bands  of  ^^Os:  Line  Posi¬ 
tions  and  Intensities,”  J.  Mol.  Spectrosc,  124, 209-217  (1987). 

16.  J.-M.  Flaud,  C.  Camy-Peyret,  V.  Malathy  Devi,  C.  P.  Rinsland, 
and  M.  A.  H.  Smith,  “The  vi  and  ^3  Bands  of  ^^O^^O^^O:  Line 
Positions  and  Intensities,”  J.  Mol.  Spectrosc.  118,  334-344 
(1986). 

17.  C.  Camy-Peyret,  J.-M.  Flaud,  A.  Perrin,  V.  Malathy  Devi,  C.  P. 
Rinsland,  and  M.  A.  H.  Smith,  “The  Hybrid-Type  Bands  vi  and 
Vs  of  ^^O^^O^^O:  Line  Positions  and  Intensities,”  J.  Mol.  Spec¬ 
trosc.  118, 345-354  (1986). 

18.  V.  Malathy  Devi,  J.-M.  Flaud,  C.  Camy-Peyret,  C.  P.  Rinsland, 
and  M.  A.  H.  Smith,  “Line  Positions  and  Intensities  for  the  + 
V2  and  Vo  +  ^3  Bands  of  ^^0$,”  J.  Mol,  Spectrosc.  125,  174-183 
(1987). 

19.  H.  M.  Pickett  et  aL,  “The  Vibrational  and  Rotational  Spectra  of 
Ozone  for  the  (0,1,0)  and  (0,2,0)  States,”  J.  Mol.  Spectrosc.  128, 
151-171  (1988). 

20.  J.-M.  Flaud,  C.  Camy-Peyret,  C.  P.  Rinsland,  M.  A.  H.  Smith, 
and  V.  Malathy  Devi,  “Line  Parameters  for  ^^Os  Bands  in  the  7- 
fim  Region,”  J.  Mol.  Spectrosc.  134, 106-112  (1989). 

21.  C.  P.  Rinsland,  M.  A.  H.  Smith,  J.-M.  Flaud,  C.  Camy-Peyret, 
and  V.  Malathy  Devi,  “Line  Positions  and  Intensities  of  the  2^3, 
vi  +  vs,  and  2vi  Bands  of  J.  Mol.  Spectrosc.  130,  204-212 
(1988). 

22.  C.  P.  Rinsland,  M.  A.  H.  Smith,  V.  Malathy  Devi,  J.-M.  Flaud, 
and  C.  Camy-Peyret,  “The  2v2  +  vs  and  2*^2  +  I'l  Bands  of  at 
4.1  fim:  Line  Positions  and  Intensities,”  J.  Mol.  Spectrosc.  139, 
343-352  (1990). 

23.  M.  A.  H.  Smith,  C.  P.  Rinsland,  V.  Malathy  Devi,  J.-M.  Flaud, 
and  C.  Camy-Peyret,  “The  3.6  i^m  Region  of  Ozone:  Line  Posi¬ 
tions  and  Intensities,”  J.  Mol.  Spectrosc.  139, 171-181  (1989). 

24.  M.  A.  H.  Smith,  C.  P.  Rinsland,  V.  Malathy  Devi,  D.  C.  Benner, 
and  K.  B.  Thakur,  “Measurements  of  Air-Broadened  and  Nitro¬ 
gen-Broadened  Half-Widths  and  Shifts  of  Ozone  Lines  near  9 
Mm,”  J.  Opt.  Soc.  Am.  B  5, 585-592  (1988). 

25.  M.  A.  H.  Smith,  C.  P,  Rinsland,  and  V.  Malathy  Devi,  “Measure¬ 
ments  of  Self-Broadening  of  Ozone,”  in  Technical  Digest,  For- 
ty-Fourth  Symposium  on  Molecular  Spectroscopy,  Ohio  State 
U.,  Columbus  (1989),  paper  TEll. 

26.  R.  R.  Gamache  and  L.  S.  Rothman,  “Theoretical  N2-Broadened 
Halfwidths  of  ^^Os,”  Appl.  Opt.  24, 1651-1656  (1985). 

27.  A.  Goldman  et  al,  “New  Spectral  Features  of  Stratospheric 
Trace  Gases  Identified  from  High-Resolution  Infrared  Balloon- 
Borne  and  Laboratory  Spectra,”  J.  Geophys.  Res.  94,  14,945- 
14,955  (1989). 

28.  F.  J.  Murcray,  J.  J.  Rosters,  R.  D.  Blatherwick,  J.  Olson,  and  D. 
G.  Murcray,  “High  Resolution  Solar  Spectrometer  System  for 
Measuring  Atmospheric  Constituents,”  Appl.  Opt.  29,  1520- 
1525  (1990). 

29.  C.  P.  Rinsland,  A.  Goldman,  F.  J.  Murcray,  F,  H.  Murcray,  and 
D.  G.  Murcray,  “Infrared  Measurements  of  Atmospheric  Gases 
Above  Mauna  Loa,  Hawaii,  in  February  1987,”  J.  Geophys.  Res. 
93, 12,607-12,626  (1988)! 

30.  M.  R.  Gunson  et  al,  “Measurements  of  CH4,  N2O,  CO,  H2O  and 
O3  in  the  Middle  Atmosphere  by  the  ATMOS  Experiment  on 
Spacelab  3,”  J.  Geophys.  Res.,  in  press. 


1  September  1990  /  Vol.  29.  No.  25  /  APPLIED  OPTICS  3671 


72 


JOURNAL  OF  GEOPHYSICAL  RESEARCH,  VOL.  94,  NO.  DI5,  PAGES  18,337-18,340,  DECEMBER  20,  1989 

APPENDIX  10 

Infrared  Emission  Measurements  of  Morning  Stratospheric  N2O5 
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Infrared  emission  spectra  obtained  during  a  balloon  flight  of  the  Air  Force  Geophysics  Laboratory 
SCRIBE  system  by  the  University  of  Denver  are  used  to  measure  stratospheric  N2O5  after  sunrise 
over  New  Mexico  (latitude  33®N).  This  is  the  first  daytime  measurement  of  N205‘  Comparisons  with 
photochemical  modeling  show  consistency  between  the  obser\^ed  and  predicted  decline  of  N2G5 
during  the  morning  hours. 


Introduction 

N2O5  is  an  important  stratospheric  reservoir  species  for 
NO.”,  produced  during  the  night  by  the  reactions 


NO2  +  O3  =>  NO3  +  O2  (1) 

NO3  +  NO2  +  M  ^  N2O5  +  M  (2) 

After  sunrise,  N2O5  gradually  photolyzes  back  into  N02^ 
N2O5  hv  ^  NO2  +  NO3  (3) 

N2O5  hv  2NO2  O  (4) 


The  diurnal  variation  of  stratospheric  N2O5,  as  computed  in 
photochemical  models  [Brassettr  and  Solomon,  1984],  is  a 
periodic  function;  it  decays  due  to  photolysis  after  sunrise 
and  increases  through  reaction  (2)  at  night. 

Nighttime  N2O5  was  tentatively  observed  by  Roscoe 
[1982]  and  Evans  [1986].  A  definitive  observation  of  N2O5 
was  made  by  Toon  et  al,  [1986]  and  Toon  [1987],  who 
observed  the  V]2  band  (^1230-1260  cm  ^  and  the  P]  and  v\\ 
bands  (~  1680-1 770  cm"^)  in  transmission  spectra  obtained 
with  the  ATMOS  instrument  at  sunrise.  Recently,  Kimde  et 
al.  [1988]  have  measured  nighttime  N2O5  from  infrared 
emission  spectra  in  the  P]2  region. 

In  the  present  work,  infrared  emission  data  are  analyzed 
to  demonstrate  the  first  measurement  of  daytime  N2O5,  and 
the  results  are  compared  with  the  predictions  of  a  one¬ 
dimensional  photochemical  model. 

Observations 

The  data  were  obtained  during  a  balloon  flight  of  the  Air 
Force  Geophysics  Laboratory  (AFGL)  SCRIBE  system 
(Stratospheric  Cryogenic  Interferometer  Balloon  Experi¬ 
ment)  by  the  University  of  Denver  from  Roswell,  New 
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Mexico,  on  July  5,  1984.  The  instrumentation  consisted  of  a 
LN2-coo1ed  Michelson  interferometer  system  employing 
cat’s  eye  optics  and  having  a  maximum  path  difference  of  ~9 
cm.  The  time  for  a  single  scan  was  ~30  s,  including  a  5-s 
fly-back.  The  field  of  view  was  -0.8°.  Data  for  the  radiomet¬ 
ric  calibration  of  the  emission  spectra  were  recorded  during 
the  flight  by  occasionally  moving  an  on-board  black  body 
into  the  instrument  field  of  view.  For  a  detailed  discussion  of 
the  SCRIBE  instrumentation,  see  Murcray  et  al.  [1984]. 

The  balloon  was  launched  at  0610  MDT  (Mountain  Day¬ 
light  Time),  and  reached  a  float  altitude  of  30.6  km  at  0740 
MDT.  Data  were  recorded  at  a  number  of  zenith  angles.  Of 
these,  the  scans  at  90.6“  and  91.1“  were  selected  as  most 
suitable  for  the  detection  and  measurement  of  daytime 
N2OJ.  The  90.6“  data  displayed  here  are  actually  the  average 
of  14  individual  scans  selected  from  those  recorded  between 
0751  and  0759  MDT,  while  the  91.1“  data  are  the  average  of 
12  scans  recorded  during  the  interval  0916-0923.5  MDT. 
Sunrise  occurred  at  the  balloon  float  altitude  at  approxi¬ 
mately  0530  MDT,  and  at  ground  level  at  the  local  time  of 
0456  MST  (Mountain  Standard  Time). 

Analysis  and  Results 

The  data  were  analyzed  by  comparison  with  synthetic 
spectra  generated  using  a  line-by-line  computer  program 
developed  at  the  University  of  Denver.  The  program,  which 
takes  refractive  effects  into  account  through  a  ray-tracing 
routine,  calculates  transmittance  or  emission  over  a  user- 
selected  number  of  atmospheric  layers  at  a  net  interval  of 
0.001  cm"'.  The  1986  edition  of  the  HITRAN  data  base 
[Rothman  et  al.,  1987]  is  used  as  input  to  the  program.  Line 
parameters  for  N2O5  do  not  exist,  but  the  absorption  coef¬ 
ficients  needed  to  compute  the  N2O5  contribution  to  the 
atmospheric  emission  are  generated  by  the  program  using 
the  absorption  cross  sections  for  N2O5  which  are  also 
included  in  the  HITRAN  data  ba.se.  Finally,  the  computed 
radiances  are  degraded  in  resolution  by  convolution  with  an 
appropriate  instrument  line  shape  function. 
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Fig.  1.  Comparisons  of  observed  (dotted  curves)  and  calculated 
(solid  cur\'es)  atmospheric  emission  spectra  in  the  1220-1320  cm~* 
region.  The  spectra  were  obtained  during  a  balloon  flight  made  on 
July  5,  1984,  by  the  University  of  Denver  with  the  AFGL  SCRIBE 
system  from  balloon  float  altitude  of  30.6  km.  Zenith  angles  and 
tangent  altitudes  are  indicated  in  the  figure.  The  top  two  frames 
include  spectral  calculations  while  for  the  bottom  frame 

N2O5  was  excluded  from  the  calculation. 


Residual  background  radiance  was  present  at  both  obser¬ 
vation  angles  (0.06  and  0.08  ftW  cm~‘  sr”Vcm”*  at  90.6° 
and  91. r,  respectively).  For  the  present  analysis  the  data 
were  adjusted  by  applying  a  wavelength  dependent  correc¬ 
tion,  described  by  a  linear  function  such  that  theory  and 
observation  agreed  in  the  ‘‘mini  windows"  centered  at 
1227.0  and  1267.07  cm~^  and  at  the  N2O-CH4  blend  at 
1268.31  cm“*.  These  three  positions  are  away  from  signifi¬ 
cant  N2O5  absorption.  The  peak  N2O5  cross  section  is  1.90 
X  10'*^  cm^  at  1246  cm (with  no  temperature  correction), 
whereas  the  cross  sections  for  the  three  calibration  points 
are  <1.5  x  lO'^^cm-.  Validity  of  the  adjustment  is  apparent 
in  the  reasonable  agreement  to  the  CH4  spectrum  longward 
of  1268  cm^*  and  insensitivity  of  the  N2O5  quantification  to 
the  background  radiance.  Subsequently,  a  least  squares 
fitting  for  N2O5  amounts  over  the  1240-1247  cm"^  interval 
(with  N2O  and  CH4  fixed)  was  made.  Recent  measurements 
of  N2O5  cross  sections  [Cantrell  et  al.,  1988/?]  show  1.79  x 
10“'*  and  1.82  X  10“'*  cm- at  1246  cm“'  for  298  and  233K, 
respectively.  Due  to  the  small  dilTerences  in  comparison  to 
the  (room  temperature)  HITRAN  value  and  the  small  tem¬ 
perature  dependence,  we  did  not  adjust  the  HITRAN  values. 

The  pressure-temperature  profile  used  in  these  simulations 
was  taken  from  a  radiosonde  ascent  from  White  Sands 
Missile  Range  at  0800  MDT  on  the  .same  day  as  the  balloon 
flight.  Mixing  ratio  profiles  for  the  various  constituents 
(except  N2O5)  are  from  Smith  [1982]. 

Figures  ]a  and  \b  show  a  comparison  between  the  calcu¬ 
lated  (solid  curve)  and  observed  (dotted  curve)  emission  for 74 
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Fig.  2.  Comparison  of  observed  (dotted  curves)  and  calculated 
(solid  cun'es)  atmospheric  emission  spectra  in  the  N2O5  1220-1270 
cm“^  region.  Figures  la  and  lb  are  expanded  versions  of  Figures 
and  Ic. 


zenith  angles  of  90,6°  and  91.1°,  respectively.  (Tangent 
heights  are  30.3  km  and  29,4  km.)  The  data  have  been 
smoothed  by  a  triangular  filter  to  0.25  cm“*  resolution.  The 
region  is  dominated  by  sharp  emission  lines  of  CH4,  N2O, 
H2O,  and,  to  a  considerably  lesser  extent,  CO2.  At  lower 
zenith  angles  there  are  also  contributions  from  HNO3  above 
—  1280  cm“^  as  well  as  from  CF4  (—1283  cm™^)  and 
CIONO2  (—1292  cm“*),  but  these  molecules  do  not  contrib¬ 
ute  appreciably  at  the  zenith  angles  shown  here.  The  broad, 
smooth  emission  feature  in  Figure  ]a  and  lb  from  1230  to 
1260  cm""^  is  due  to  stratospheric  N2O5.  This  can  be  seen 
from  Figure  Ic,  which  is  the  same  as  Figure  la^  except  that 
no  N2O5  has  been  included  in  the  calculation.  Figure  2 
shows  an  expanded  plot  of  the  N2O5  1220-1270  cm  region 
from  Figures  1^2  and  Ic.  Figure  3  presents  the  observed 


1220  1230  1240  1250  1260  i270 


WAVENUMBER  (cm-’) 

Fig.  3.  Observed  emission  spectra  (dotted  curves;  ana  theoretical 
calculation  (solid  curve)  only  including  the  N2O5  contribution. 
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Fig.  4.  One-dimensional  photochemical-transport  model  (see 
text)  mixing  ratio  profiles  (solid  curves)  of  N2O5  for  predawn  (0440 
MST)  and  morning  conditions  used  for  comparison  with  the  exper¬ 
iment.  The  triangles  show  the  measured  N2O5  mixing  ratios  1.26 
and  1.04  ppb,  at  the  tangent  altitudes  30.3  and  29.4  km. 


radiance  at  90.6"  and  a  theoretical  calculation  which  only 
includes  the  contribution  due  to  N2O5. 

The  synthetic  spectra  of  Figures  la  and  \b  were  computed 
using  N2O5  profiles  of  the  same  general  shape  as  the  solid 
curves  of  Figure  4,  which  are  generated  as  described  below, 
multiplied  by  a  scaling  factor.  The  triangles  in  Figure  4  show 
the  resulting  N2O5  mixing  ratios,  1.26  and  1.04  ppb,  at  the 
tangent  altitudes  of  the  two  scans,  30.3  and  29.4  km, 
respectively.  Theoretical  N2O5  mixing  ratios  are  2.10  and 
1 .65  ppb  for  the  corresponding  times  of  0655  and  0820  MST. 
The  observed  relative  decrease  in  N2O5  is  therefore  0.83, 
whereas  the  theoretical  value  is  0.79.  The  estimated  accu¬ 
racy  of  the  absolute  N2O5  measurements  is  ±45%.  This  is 
based  on  uncertainties  in  radiance  calibration  and  pointing 
(10%),  vertical  slope  of  the  theoretical  N2O5  mixing  ratios 
(10%),  spectral  line  parameters  (15%),  and  the  N2O,  CH4 
profiles  (10%).  Error  in  the  observed  relative  change  is 
estimated  as  ±  15%  (due  to  partial  cancellation  of  systematic 
errors). 

The  theoretical  N2O5  mixing  ratios  curves  were  calculated 
by  a  one-dimensional  photochemical  transport  model  [Cice¬ 
rone  et  al.,  1983].  The  time  dependent  calculation  extended 
between  10  and  80  km  altitude.  The  pressure-temperature 
profile  was  that  mentioned  above,  and  the  solar  declination 
of  -i-23"  and  the  latitude  of  +33°  matched  that  for  the  balloon 
flight. 

Reaction  rates  for  81  two-  and  three-body  gas-phase 
processes  and  25  photolysis  processes  for  38  gaseous  species 
follow  that  of  the  Jet  Propulsion  Laboratory  (JPL  1987) 
[DeMore  et  al.,  1987],  but  with  the  Cantrell  et  al.  [1988i>] 
equilibrium  constant  for  reaction  (2).  At  30  km  at  a  temper¬ 
ature  of  23 IK,  the  Cantrell  and  JPL  1987  equilibrium  con¬ 
stants  differ  by  27%.  However,  for  the  geometry  of  the  July 
5  observation,  the  rate  of  thermal  decomposition  of  NjOs  is 
slower  than  the  diurnal  average  of  N2O3  photolysis  by  a 
factor  of  78,  so  the  influence  of  the  revised  equilibrium 
constant  upon  the  N2O5  mixing  ratio  is  small.  The  total 
odd-chlorine  (Cl  +  CIO  +  HOCl  +  CIONO2  +  HCl  +  CIO2) 
mixing  ratios  are  2.17  and  2.61  ppb  at  30  and  60  km, 
respectively.  The  J  values  at  30  km  are  3.42  x  10  *  s  '  and 
5.78  X  10"^  s“'  forN205  photolysis  (reactions  (3)  and  (4)) 
for  0655  and  0820  MST,  respectively.  The  albedo  value  was- 


0.35,  and  J  is  insensitive  to  changes  in  the  albedo  (e.g.,  J 
equals  5.62  x  10"^  s"’  for  an  albedo  of  0.25  at  0820  MST). 

Since  N->05  is  but  one  member  of  the  NO,,  family  (N  + 
NO  +  N62  +  NO3  +  HNO3  +  HNO4  +  2*N205  + 
CIONO2),  the  theoretical  value  of  N2O5  is  related  to  the 
total  n6^  mixing  ratio.  For  the  theoretical  curves  of  Figure 
2,  NOy  is  equal  to  15.0  and  14.3  ppb  at  40  and  30  km, 
respectively.  Obserx'ations  of  Russell  et  al.  [1988]  for  May 
1985,  sunset,  30°N  latitude  indicate  NOy  near  16.7  and  15.2 
ppb  at  40  and  30  km,  roughly  11%  higher  than  our  one¬ 
dimensional  theoretical  predictions. 

The  model  N2O5  and  NOy  profiles  are  reasonable,  though 
they  may  not  accurately  represent  the  mixing  ratio  profiles  of 
July  5,  1984,  due  to  seasonal  variations  and  model  limita¬ 
tions.  Differences  between  one-  and  two-dimensional  models 
are  illustrated,  for  example,  in  Figure  10-58  of  World  Mete¬ 
orological  Organization  (WMO)WASA  [1986],  in  which  the 
two  modeling  techniques  at  30  km  estimate  different 
amounts  of  NO,,  and,  presumably,  different  amounts  of 
N2OJ.  Despite  such  differences,  the  present  measurement  of 
the  relative  change  of  N2O5  during  the  morning  hours  is 
consistent  with  the  temporal  change  of  the  one-dimensional 
calculation.  Our  observed  early  morning  decline  in  N2O5 
concentration  is  also  consistent  qualitatively  (see  equations 
(1)_(4))  with  measured  early  morning  increases  in  both  NO 
and  NO2  [e.g.,  Ridley  et  al,  1977;  Flaud  et  al.,  1988; 
Rinsland  et  al.,  1988]. 

The  present  morning  N2O5  results  and  the  sunrise  and 
nighttime  N2O5  observations  by  Toon  et  al.  [1986]  and 
Kunde  et  al.  ’[1988],  all  report  N2O5  mixing  ratios  in  the  1-  to 
2-ppb  range  near  30  km  altitude.  One  sees  that  our  obserx'ed 
mixing  ratios  at  0755  MDT  and  0920  MDT  are  -60%  smaller 
than  the  values  predicted  by  our  model.  Kunde  et  al.  [1988] 
showed  that  their  observed  nighttime  N2O5  is  20-30% 
smaller  in  comparison  to  a  (different)  photochemical  model. 
However,  the  relative  decrease  of  0.83  over  the  1.5  hours 
between  the  observations  is  consistent  with  that  of  0.79 
predicted  by  the  model.  The  relative  decrease  in  N2O5 
concentration  during  the  morning  can  be  approximated  by 
D  =  exp  (— /  A/),  where  J  is  the  average  N2O5  photolysis 
rate  over  the  time  interval  Ar.  With  J  =  4  x  10“^  s  '  and  At 
=  1  hour  25  min,  D  =  0.81,  in  good  agreement  with  the 
detailed  model  results  of  0.79.  D  is  not  sensitive  to  the  actual 
N2O5  amount,  and  thus  the  absolute  difference  between 
models  and  measurements  of  N2O5  is  less  important  than  the 
temporal  variations. 
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